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ABSTRACT 


A  general  repairable  system  model  is  proposed  with  a 
maintenance  schedule  such  that  the  failure  rate  is  reduced 
after  each  preventive  maintenance.  The  model  introduces 
naturally  the  concept  of  system  availabilities  as  random 
variables.  Probability  distributions  of  the  stochastic 
availabilites  are  derived.  The  stochastic  availability  is 
optimized  with  respect  to  the  duration  of  the  operating 
interval  between  maintenance  operations,  and  with  respect 
to  an  age  replacement  policy.  Analysis  and  optimization  are 
achieved  by  a  PASCAL  computer  program  which  is  developed  for 
computations  and  illustrations  of  various  failure  rate  charac¬ 
teristics  and  parameter  variations.  The  concept  is  further 
extended  to  stochastic  cost  function  or  generalized  stochastic 
avai lability. 
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CHAPTER  1 


INTRODUCTION 

§1.1  CONCEPTS  OF  SYSTEM  AVAILABILITY 

Reliability  is  a  relatively  new  field  of  engineering  (Shooman  (1968); 
Dhillon  and  Singh  (1981))  and  applied  probability  (Barlow  and  Proschan  (1965)). 
It  has  found  its  importance  in  the  planning,  design,  and  operation  of  systems 
in  recent  years.  Formally,  system  reliability  is  defined  as  the  probability 
that  a  system  operates  without  failure  during  the  time  interval  CO,t]  , 
that  is, 

P{  T >  t  }  «  1  -  F(t)  (1.1.1) 

where  T  is  the  (random)  lifetime  of  failure-free  operation  of  an  initially 
good  system  until  failure,  F(.)  ;  is  the  distribution  function  (d.f.)  of  T 
and  p{.}  denotes  the  probability  measure.  The  reliability  is  mainly 
defined  for  unrepairable  systems.  Other  definitions  can  be  found  in  various 
references  (Barlow  and  Proschan  (1965);  Kozlov  and  Ushsakov  (1970)). 

Note  that  the  field  of  reliability  engineering  concerns  the  probability  of 
non-negative  random  variables  only. 

For  repairable  and  maintained  systems,  a  measure  of  system  performance 
and  a  main  design  crierion  is  the  concept  of  system  availability.  Various 
definitions  of  availability  have  been  defined  in  the  past  three  decades 
(Osaki  and  Nakagawa  (1976);  Lie,  Hwang  and  Tillman  (1977)).  We  shall 
review  some  of  them. 

Hasford  (I960),  Barlow  and  Hunter  (1960),  and  also  Barlow  and  Proschan 
(1965)  defined  the  following  three  types  of  availability  : 
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(1)  Instantaneous  Availability  /  Pointwise  Availability 

The  instantaneous  (or  pointwise)  availability  is  the  probability  that 
the  system  is  operative  at  any  time  t  .  It  gives  a  good  mesure  for 
systems  functioning  at  any  random  time. 

(2)  Average  Uptime  Availability  /  Interval  Availability 

The  average  uptime  (or  interval)  availability  is  the  expected  ratio  of 
uptime  in  a  given  time  interval.  It  is  a  satisfactory  measure  for 
system  working  over  a  duty  cycle. 

(3)  Steady-State  Availability  /  Limiting-Interval  Availability 

The  steady-state  (or  limi ting- Interval )  availability  is  the  average 
uptime  availability  when  the  time  interval  is  very  large.  It  is 
relatively  simple  to  calculate  and  it  is  a  quite  satisfactory  measure 
for  continously  operating  systems.  (See  also  Eq.(1.1.3)). 


Kabak  (1969)  proposed  the  following  two  versions  of  availability  for 
systems  with  up  and  down  cycles  for  an  exponential  failure  time  and  constant 
repair  time  : 

(a)  Availability  for  Multipe  Cycles 

The  availability  of  n  cycles  (  n  ■  1,2,3,...  )  is  the  expected  value 
o‘f  the  proportion  of  total  uptime  in  the  n  cycles  to  the  total  elapsed 
time  in  the  n  cycles. 

(b)  Finite  Time  Availability 

The  finite  time  availability  A(T)  for  a  time  interval  CO,T)  is 
determined  by  combining  the  probability  of  n  failures  in  a  given  time 
■  interval  with  the  proportion  of  available  time  for  the  interval. 

(See  Kabak  (1969)  for  detail). 


3 


It  has  been  shown  by  Kabak  (1969)  that  for  exponential  failure  time  and 
constant  repair  time, 

(a)  as  the  number  of  cycles  increases,  the  limit  of  the  availability  for 
multiple  cycles  is  the  steady-state  availability,  and 

(b)  as  the  time  interval  increases  the • finite ■ time  availability  also 
approaches  the  steady-state  availability. 


Consider  a  repairable  system  with  up  and  down  cycles  all  of  which  are 
independent.  Let 


T  a  uptime  random  variable  in  an  up-cycle 

up 

T ,  s  downtime  random  variable  in  a  down-cycle 
down 


T ' s  are  independent  and  identically  distributed  ( i . i . d . ) 


Note  that  all  the 

up 

random  variables  and  all  the  Td0wn ' 3  a**8  also  i.i«d.  It  can  be  shown  by 
renewal  theory  for  an  alternating  renewal  process  (Parzen  (1962))  that 


lira  p{  System  is  up  at  t  ]■ 
t -to 


*(Tuo) 

EfT  1  +  e(t.  1 
l  upJ  *  down* 


MTBF. 

MTBF  +  MTTR 

where  E[.}  3  expectation 

MTBF  s  mean  time  between  failure 
MTTR  3  mean  time  to  repair. 


(1.1.2) 


(1.1.2') 


The  fraction  in  Eq.(1.1.2)  is  simpy  the  steady-state  availability  which 
will  be  denoted  by 

A® 

Thus, 


E(  uptime  3 _ 

EC  uptime  }  +  E(_  downtime  ) 


(1.1.3) 


In  particular,  for  exponentially  distributed  uptime  and  downtime,  namely, 
with  probability  density  functions  (p.d.f.) 
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fT  <t)  *  A«xp(-At)  ,  t^o 

and 

fT 

down(t)  =  ^»exp(-/tt)  ,  t^O 

then 

A®  „  1/A  =  M 

1/A  +  1/m  /“  +  A 


(1.1.4) 


(1.1.5) 


(1.1.6) 


Eq. (1.1.3)  is  widely  accepted  as  the  (simplest  and  non-random) 
definition  of  system  availability.  However  Eq. (1.1.3)  gives  only  the 
average  value  and  there  is  no  probabilistic  guarantee  that  A®  will  ever 
be  achieved.  A  partial  remedy  had  been  proposed  by  Martz  (1971)  who 
defined  the  following  random  variable  : 


A  3 


uptime 


uptime  +  downtime 


(1.1.7) 


JiE. 


T  +  T , 
up  down 


(1.1.7') 


Since  T  and  T.  are  random  variables,  A  is  indeed  a  random  variable, 
up  down 

Martz  studied  the  following  definition  of  availability  : 


Single  Cycle  Availability 

The  single  cycle  availability  Ar  is  the  value  such  that 
p{A>Ar}  -  r 
for  all  r£(0,l)  ,  or 


1 


f.  (a)  da  *  r 
A 


(1.1.8) 


(1.1.8') 


where  f. (.)  is  the  p.d.f.  of  A  . 
A 


By  specifying  r  ,  Marts* s  definition  gives  a  probabilistic  guarantee 
on  the  frequency  of  occurrence  of  the  availability  value  .  Marts  also 
derived  some  expressions  of  A  by  evaluating  f  ( . )  for  independent 

”  A 

uptime  and  downtime.  In  particular,  if  T  and  Tdown  are  exponential 
as  given  in  Eq. (1.1.4)  and  (1.1.5},  then 


A  .  - ■  (1.1.9) 

( l-r)/i  +  tX 

When  ra%  , 

K  -  A®  (1.1.10) 

'Z 


Some  complicated  expressions  for  Ar(t)  had  been  derived  by  them  for 
independent  uptime  and  downtime. 
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Marshall  and  Goldstein  (1980)  studied  a  repairable  system  and 
introduced  the  concept  of  cycle  availability  as  a  random  variable  : 

Cycle  Availability 

The  cycle  availability  A  is  a  random  variable  defined  a3 

c 


Total  uptime  for  a  complete  cycle 

.Total  time  for.  .Total  downtime  for. 
la  complete  cycle'  +  '  a  complete  cycle 


(1.1.13) 


where  a  complete  cycle  of  the  system  begins  with  a  new  operative  system 
and  ends  until  failure,  followed  by  a  renewal  or  complete  repair. 


Goldstein  considered  system  designs  and  studied  the  following 
probabilistic  inequality 

P  {ac?  1— 1|  ^  1-5  (1.1.14) 

where  0<t<l  and  0< $<  1  .  He  derived  some  expressions  for  repairable 
systems  with  independent  exponential  lifetime  and  exponential  repair  time, 
so  that  the  technique  of  Laplace  transform  was  applicable  to  give  an  Erlang 
distribution  (Parzen  (1962)).  However  no  numerical  illustrations  were 
given.  The  generalization  of  independent  exponential  distribution  to  more 
general  distributions,  such  as  Weibull  distribution  (see  Appendix  1),  is 
not  ‘straight  forward,  because  the  Laplace  transform  technique  does  not  give 
a  closed  form  solution. 


In  this  research,  various  concepts  of  system  availability  are  studied 
as  random  variables,  and  they  will  be  called  stochastic  availabilities,  in 
order  to  emphasize  their  random  nature.  In  particular,  Goldstein's  model 
and  approach  are  extended  to  a  general  repairable  system  with  maintenance 
schedule  such  that  the  failure  rate  is  both  age  and  maintenance  dependent. 
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§  1.2  MAINTENANCE  AND  FAILURE  RATE 


As  we  have  seen  so  far,  the  various  definitions  and  expressions  of  system 
availability  have  no  explicit  dependence  on  the  maintenance  or  repairs. 

Most  authors  have  treated  failure  and  repair  as  independent  random  variables, 
neither  of  which  depends  on  the  maintenance.  In  other  words,  the  system 
failure  rate  is  not  affected  by  any  preventive  maintenance.  The  use  of 
failure  rate  seems  to  be  a  natural  way  of  system  failure  analysis  because  a 
couplete  knowledge  is  obtained  once  we  know  or  specify  the  failure  rate 
(see  Theorem  2.2.1). 

The  purpose  of  preventive  maintenance  for  repairable  system  is  most 
likely  to  improve1  the  system  if  possible.  However  the  failure  rate  may  be 
disturbed  by  the  number  of  maintenance  operations.  Various  replacement 
and  maintenance  policies  had  been  proposed  (Barlow  and  P roach an  (1965); 
Jorgenson,  McCall  and  Radner  (1967);  and  Pierskalla  and  Voelker  (1976)). 
Basically,  an  optimum  replacement  policy  is  to  choose  a  set  of  time  intervals 
between  two  maintenance  periods  at  which  replacements  are  to  take  place 
such  as  to  minimize  an  expected  cost  during  a  given  finite  or  infinite  time. 
Note  that  the  cost  is  a  deterministic  quantity.  Recently,  Nguyen  and  Murthy 
(1981)  have  extended  Barlow  and  Hunter's  (1960)  two  replacement  policies  to 
a  case  where  the  failure  rate  is  an  increasing  function  with  the  number  of 
repairs.  A  further  generalization  of  the  optimization  problem  to  failure 
rate  dependence  on  the  number  of  previous  repairs  and  on  the  times  when  they 
took  place  was  proposed  by  Shaw,  Sbrahimian  and  Chan  (1981). 

As  with  availability,  the  cost  to  be  minimized  can  be  generalized  to 
a  random  variable  approach  by  modifying  the  definitions  of  stochastic 
availability  to  be  associated  with  costs,  so  that  minimizing  the  probability 
of  the  so  called  stochastic  cost  is  equivalent  to  maximizing  that  of  the 


stochastic  availability  (see  Chapter  5) 


§1.3  SUMMARY 


In  this  research  various  concepts  of  stochastic  availability  for  a 
general  repairable  system  model  with  a  maintenance  schedule  and  an  age  and 
maintenance  dependent  failure  rate  are  studied  and  are  applied  to  system 
design.  In  Chapter  2,  the  general  system  model  is  proposed.  The  system 
failure  rate  is  reduced  and  altered  after  each  preventive  maintenance. 

Two  types  of  failure  reduction  are  considered.  The  lifetime  distributions 
for both  failure  reduction  types  are  derived.  In  Chapter  3,  various 
stochastic  availabilities  are  defined.  The  probability  distributions  of 
the  stochastic  availabilities  an  derived.  In  Chapter  4,  optimum  system 
designs  are  considered  using  the  concept  of  maximizing  the  probability  of 
stochastic  availability  and  probabilistic  inequality.  In  particular, 
optimum  operative-cycle  time  and  age  replacement  time  are  determined  for  a 
periodic  operative-maintenance  policy.  Numerical  examples  illustrate  the 
failure  rate  characteristics  and  variations  of  system  parameters  and  age 
replacement  time.  A  PASCAL  computer  program  is  developed  for  the  computation 
of  the  probability  distributions  of  stochastic  availabilities.  In  Chapter  5f 
the  concept  is  further  extended  to  the  so  called  generalized  stochastic 
availability  and  the  equivalent  stochastic  cost  function  in  probabilistic 
optimization.  In  Chapter  6,  the  work  is  concluded  with  a  discussion  on 
possible  extensions.  Appendix  1  gives  a  summary  of  Veibull  distribution. 
Appendix  2  is  a  listing  of  the  PASCAL  computer  program. 


CHAPTER  2 


SYSTEM  MODEL 


§2.1  A  REPAIRABLE  SYSTEM  WITH  MAINTENANCE  SCHEDULE 


We  shall  study  a  repairable  syttem  with  a  maintenance  schedule  as 
shown  in  Fig.  2.1.1. 


Maintenance 

State 


Fig.  2.1.1  State  transition  diagram  of  a  repairable  system 
with  a  maintenance  schedule. 
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The  system  starts  Its  mission  when  it  is  new  and  it  follows  an  operative 
and  maintenance  schedule  before  it  fails.  Thus  the  system  is  represented  by 
an  UP  operative  state  @  of  normal  operation,  a  DOWN  maintenance  state  (9) 
of  preventive  maintenance.  By  preventive  maintenance  we  mean  that  the  system 
is  not  replaced  by  a  new  one  completely,  but  simply  a  replacement  of  parts, 
minors  repairs  or  even  just  inspections.  A  DOWN  failure  state  ©  accounts 
for  the  system  failure  or  breakdown.  A  new  system  is  set  to  operate  for  a 
time  period  called  the  operative-cycle  time  (u)  and  the  system  is  said  to  be 
in  the  operat i ve-cy c 1 e  (or  simply,  the  0-cycle) .  Then  it  is  brought  to 
preventive  maintenance  for  a  duration  called  the  maintenance-cycle  time  (d) 
and  the  system  is  in  the  maintenance-cycle  (the  M-cycle) .  After  the 
maintenance  operation  the  system  is  back  to  normal  operation.  The  transition 

is  called  an  operative-maintenance  cycle  (an  0-M  cycle) .  The  system  may 
fail  during  the  0-cycle  or  during  the  M-cycle.  When  the  system  breaks  down, 
it  is  in  the  failure  state  ©  ,  and  it  will  only  resume  its  normal  operation 

as  a  new  system  after  a  complete  repair,  renewal  or  replacement  of  the  entire 
system.  The  time  spent  for  such  a  renewal  is  assumed  to  be  a  random  variable 
called  the  repair  time  Tp  .  If  the  transition 

©—  <g) 

4 

is  completed,  the  system-  is  said  to  have  gone  through  a  complete  cycle  and 
the  entire  system  is  considered  to  be  a  new  system  again. 

A  schedule  of  the  0-M  cycles,  operative-cycle  time  and  Maintenance-cycle 
time  is  shown  in  Fig.  2.1.2. 
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Tine  : 

State  : 

Cycle  : 

Fig.  2.1.2  System  operative-maintenance  schedule  before  failure. 


We  have  used  the  following  notations  in  Fig.  2.1.2  : 


v° 

0-Cycles  : 

CW  . 

^l’V . ^n-l’V 

M-Cycles  : 

• 

Cc1»t1)  * 

CW . ’  V 

0-M  Cycles  : 

*vv  * 

*W . ^n-l’V 

Operative-cycle  tine  : 

WS  • 

...  t  t  .  ,  ... 

n  n  n— i 

Maintenance-cycle  time  : 

W*i  . 

*'•  '  W*n  '  *•* 
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In  particular,  we  shall  mainly  interest  in  the  periodic 
operative-maintenance  schedule  (see  also  §  4.1),  namely, 


U,  ■  U_  m  ...  »  U  a. ..a  u 

12  n 

d.  ■  d_  ■  • . .  ■  d  ■*..  a  d 
12  n 


t 


n 


n.t 


1 


t^  ■  u  +  d 


>  (2.1.1) 


NOTE  :  All  time  variables  (deterministic  or  random)  are  assumed 

non-negative.  Therefore,  all  probability  density  functions 
and  failure  rate  functions  are  defined  to  be  zero  for  any 
negative  time.  We  shall  not  specify  this  assumption  explicitly 
in  all  expressions. 


13 


§2.2  FAILURE  RATE  WITH  AGE  AND  MAINTENANCE  DEPENDENCE 


The  lifetime  of  the  repairable  system  before  failure  is  a  random  variable 
T  with  a  distribution  function  (d.f. )  F(.)  known  as  the  failure  distribution 

function  (Barlow  and  Proschan  (1965)).  Usually,  F(.)  is  assumed  to  be 
absolutely  continuous  so  that  there  exists  a  probability  density  function 
(p.d.f.)  given  by 

««  -  ^ 

or,  ft 

F(t)  S  J0  f(x)  dx 


Note  that 

f(t)  -  0 

for 

t<0 

and 

F(t)  »  0 

for 

t<0 

system  hazard  rate  function 


(2.2.3) 


(2.2.4) 

(2.2.5) 


The  system  failure  rate  function  (f.r. )  or  the 
r( . )  is  defined  by 


r(t)  a 


f(t) 


1  -  F(t) 

dR(t) 

dt 


where 


R(t) 

R(t)  a  l-  F(t) 


F(t) 


(2.2.1) 

(2.2.2) 


is  called  the  system  reliability  function  (r.f.). 
From  Eq. (2.2.4),  we  have  • 

+  r(t)R( t)  -  0 


1 

i 


Solving, 


R(t)  *  1  -  F(t)  -  exp |  -  p  r(x)  dx 

R(0)  -  1 

F(0)  -  0 

Eq. (2.2.6)  and  (2.2.3)  give 

f(t)  -  r(t)R(t) 

or 

f(t)  m  r(t).exp  p  r(x)  dx  | 

Eq. (2.2.7')  is  of  fundamental  importance. 

In  order  to  be  a  valid  failure  rate  function,  r(t)  must  satisfy  the 
following  conditions  : 

r(t)  «  0  ,  t<0 

r(t)  £  0  ,  t£0 

ao 

r(x)dx  "  oo 

0 

From  Eq. (2.2.2),  (2.2.3)  and  (2.2.7*),  we  see  that  F(.),  f(.)p  and  r(.) 
are  equally  suitable  for  describing  the  failure  distribution.  In  particular, 
if  we  specify  the  failure  rate,  then  the  following  theorem  is  obvious  : 

THEOREM  2.2.1 

A  failure  rate  function  r(.)  satisfying  (2.2.8) 
uniquely  determines  the  failure  distribution  of  the 
system. 


(2.2.7) 

(2.2.7*) 


1L. 


Since  there  are  infinitely  many  failure  laws,  we  shall  choose  a  widely 
useful  one,  namely,  the  Weibull  type  failure  rate  (Appendix  1) ,  that  is, 
the  system  lifetime  T  has  a  Weibull  distribution  : 


P.d.f. 

o 

f(t)  - 

jf^t*-1  exp(-J?t-) 

(2.2.9) 

D.f  . 

S(t)  - 

1  -  exp(-jft“) 

(2.2.10) 

F.r. 

2(t)  - 

(2.2.11) 

R.f. 

fi(t)  = 

exp(-Aa  tf* ) 

(2.2.12) 

In  particular, 

(1)  0(  -  l  (exponential  case),  constant  failure  rate; 

(2)  a  m  2  (Rayleigh  case),  linear  failure  rate. 

Therefore,  the  failure  rate  is  age  dependent  whenever  o(  jt  1  . 

We  shall  study  the  Rayleigh  case  (  01  «  2) ,  or  the  linear  failure  rate, 
because  all  other  cases  can  be  analyzed  in  the  same  way. 

In  order  to  allow  the  failure  rate  r(t)  to  be  dependent  on  the 
maintenance  operations,  the  failure  rate  is  assumed  to  change  after  each 
maintenance-cycle.  Thus  the  failure  rate  is  piecewisely  defined  on  each 
0— M  cycle  in  the  following  manner  : 

,  t  € (tQ,ti) 

,  t€(tx,t2) 

(2.2.13) 

,  tec*  ,,t ) 

’  v  n-1  n 


where 

rx(t),  r2(t),  ...  ,  rR(t),  ... 

is  a  sequence  of  Weibull  type  failure  rates  subjected  to  (2.2.3)  ,  namely, 


ICi 


f  r2(t)^0  ,  ...  ,  rn(t)^0  , 


*1 

’*2 

r^xjdx  + 

r2(x)dx  +  ...  + 

*0  J 

r(x)dx 

n 


+  . .  • 


QO 


The  following  form  of  the  failure  rate  is  proposed  : 


For  n  »  1,2,3,  ... 


rn(t)  -  ?(t;>n,«)  -  A(t)  ,  t^CVi’V 


where 

0  ^  A  *  ^2^  ^3^  •  •  •  $  An  ^  ^  CO 

^(t)  -  »  t^^tn^l,tn) 

0  -  &i  <  Aa  ^  —  4  An^  - 

A„  -  A,,.!  *  Wl'Vl1 
81  ’  0  ’  *n€  to,1J  ’  °  "  2’3’-” 


(2.2.14) 


(2.2.15) 

(2.2.16) 

(2.2.17) 

(2.2.18) 

(2.2.19) 

(2.2.20) 
(2.2.21) 
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The  failure  rate  r(t)  so  defined  has  discontinuities  of  "jump"  type  at 
a  countable  set  of  points 


*0  '  \ . tn  *  '** 

and  A (t)  also  has  discontinuities  at  the  same  points  representing  a 
reduction  failure  rate.  Obviously,  if  there  is  no  failure  reduction,  then 


^(t)  -  0  — »  r(t)  »  r(t;A,«) 


Thus  the  benefit  of  the  maintenance  schedule  improves  an  aging  system  in  the 
sense  that  the  system  failure  rate  i3  reduced  after  each  0-M  cycle  with  a 
tradeoff  of  shortening  the  mean  lifetime  because  the  X  ' s  which  is  a  measure 
of  mean  lifetime  (Appendix  1,  Eq.(A1.4)),  increase  with  the  number  of 
O-M  cycles.  For  convenience,  we  define  the  following  : 


r(t;A,«) 


regular  failure  rate  (3  ?(t)  if  /  and  Of  are  constants) 


r(t) 

A(t) 


age  and  maintenance  dependent  failure  rate 
failure  reduction  function 
failure  reduction  factor 


An  a  reduction  jump-down 


As  an  illustration  we  depict  the  shapes  of  r(t),  ^(t),  and  r(t)  in 
Fig.  2.2.1  for  a  Weibull  type  linear  failure  rate  (ec»  2). 

When  £^(t)  is  a  step-like  function,  the  r(t)  has  a  constant  slope 
in.  each  interval.  When  ^(t)  is  a  piecewise  linear  function,  the  r(t) 
ha3  a  varying  slope  in  each  interval.  This  corresponds  to  a  maintenance- 
dependent  failure  rate  parameter  sequence  f  }  . 


la 


|  2.3  EAILURE  REDUCTION  CRITERIA 


For  a  real,  usually  large,  system,  it  is  more  efficient  to  perform 
preventive  maintenance  on  a  periodic  schedule  in  order  to  reduce  management 
cost  etc.  It  is  therefore  reasonable  to  assume  that  the  operative— cycle 
times  and  maintenance-cycle  times  are  fixed  (Eq. (2.1.1) ) , 


U1  *  u2 
dl-d2 


...  ad  3  ...  3  d 
n 

t^  a  ,  &  a  1,2,3,... 


Furthermore,  we  may  assume  that 


and 


A,  "  =  ...  *  A  =>.. .»  A 

i  £.  n 


*  S2  3  •••  ■  gn  a  ...  s  g  £  [OflJ 


f  (2.3.1) 


(2.3.2) 

(2.3.3) 


Two  special  types  of  failure  reduction  are  defined  in  the  following  : 


TYPE  1  FIXED  FAILURE  REDUCTION 


After  each  M-cycle  the  failure  rate  is  reduced  in  a  way  that  all 
jump-downs  are  the  same  in  each  Q-M  cycle,  or  more  precisely, 


-  0 

A2  ■  s-VV  £  A 
A3  -  A2  +  A  -  2A 

A„  -  An_i  +  A  -  (n-D.A 


\ 


>  (2.3.4) 


where 


g£Co,lJ 


J 
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Thus, 


X  -  2 

— ►  A 

2  2*\tl 

(2.3.5) 

The  failure  rate  is  given  by 

’  2A*t 

,  tSCtQ,^) 

2A*t  -  £ 

,  t  €  ltvt2) 

r(t)  -  * 

2A*t  - 

• 

,  t6£t2,t3) 

(2.3.6) 

• 

2  A*t  -  (n-l)A 

e 

• 

• 

•  tftWV 

The  shape  of 

r(t)  is  sketched 

in  Fig.  2.3.1.  The  r(t) 

given  by 

Eq. (2.3.6)  will  be  called  the  fixed  reductuion  piecevrLsely  linear  failure  rate. 


The  case  when  the  failure  rate  parameter  A  is  not  constant  is  also  of 
interest.  We  shall  assume  that  A  is  actually  maintenance-dependent, 
namely,  A  varies  in  the  form  of  an  increasing  sequence  as  in  (2.2.17). 

This  is  the  case  when  the  mean  lifetime  of  the  repaired  system  decreases 
with  the  number  of  preventive  maintenance  carried  out.  (From  Eq.(Al,4), 
the  mean  lifetime  is  inversely  proportional  to  the  A'8»)  We  shall 
illustrate  this  case  in  Fig.  4.4.8  where  the  Vs  are  defined  by  Sq.  (4.4.4). 


Fig.  2.8.1  Fixed  reduction  linear  failure  rate 
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TYPE  2  PROPORTIONAL  FAILURE  REDUCTION 

After  each  M-cycle  the  failure  rate  is  reduced  such  that  each  "jump-down" 
is  proportional  to  the  value  of  the  failure  rate  at  that  instant,  namely. 


-  o 

A2  •  4-W 

■  g.r2(t2)  +  g.r^tj.) 

An  *  «•**„(  tQ)  +  •••  + 

e 

where 

g6Co.il 


(2.3.7) 


Thus,  for  OL  m  2  , 

r^t)  »  2>xt 

r2(t)  «  2jJt  -  2A*gtx 

r3(t)  *  2  Xt  -  EA^gttg+d-gJt^) 


(2.3.8) 


The  shape  of  r(t)  is  sketched  in  Fig.  2.3.2.  The  r(t)  given  by 
Eq. (2.3.8)  will  be  called  the  proportional  reduction  piecewisely  linear 
failure  rate. 


ito 


::a 


2.4  RELIABILITY  FUNCTIONS 


Rocall  that  the  system  reliability  function  U(t)  is  related  to  the 
failure  rate  function  r(t)  by  Eq. (2.2.6). 

Define 

r(t)  »  A*aft*_1  a  regular  failure  rate  (2.4.1) 

o  ft 

R(t)  a  exp(-  J0  r(x)dx)  s  regular  reliability  function  (2.4.2) 


For  an  age  and  maintenance  dependent  failure  rate  r(t)  defined  by 
(2.2.13)  -  (2.2.21),  we  define 


R(t)  3  ex p(-  J**  r(x)  dx  ) 


and  for  t  >t  , 
*  n— l 


n-1 

rr  r  (t ) 

j»i j  j 


exp( 


-n 


r(x)dx  ) 


n-1 


where 


R.(t)  5  exp(  -  f  r(x)dx  )  ,  t  €  [t  ,,t.)  ,  j»l,2,, 
j  j-1  J 

If  all  's  in  (2.2.17)  are  constant,  then 

rn(t)  -  ?(t)  -  ^ 

which  ,  together  with  Eq. (2.4.5),  give 


(2.4.3) 


R(t)  a  exp  (  -  j^rMdx  -  j^2r(x)dx  -  ...  -  j  '  r(x)dx  )  (2.4.4) 

i  n-1 


(2.4.4*) 


(2.4.5) 


n^(t)  a  exp(  - 

(J  r.(x)dx  ) 

}tj-l  J 

a  exp(  - 

^  r(x)dx  ).exp(  [t-t^^Aj) 

j  * “1 

(2.4.6) 

Therefore,  Eq. (2.4.4') 

bocomes 

o 

R(t)  a  R(t).exp 

[f  ‘VVl^d  *  ^V^An] 

(2.4.7) 

4-1 

t 


(1)  FIXED  REDUCTION  RELIABILITY  FUNCTION 


I 


r 


for  fixed  reduction,  (2.3.4)  and  (2.4.7)  gives 

R(t)  -  R(t).o  *{[£  (j-l)(tj-tj_1)  +  (n-l)  ( ‘t~tn_1)  ]a]  «..*•« 

j-1 


For  constant  operative-cycle  and  maintenance-cycle  times  (periodic 
operative-maintenance  schedule), 


t  *  n.t, 
n  1 

(2.4.9) 

we  have. 

0  r  i 

R(t) 

»  R(t).exp£  (n-l)(t  -  J^nt^)^  J 

(2.4.10) 

<*  «  2  — ^ 

R(t)  ■  exp(  -A?t2  +  2(n-l)(t-/^nt1)g  A*t.) 

(2.4.11) 

(2)  PROPORTIONAL  REDUCTION  RELIABILITY  FUNCTION 

For  proportional  reduction,  with  'r(t)  given  by  Eq. (2.3.9),  we  have. 


R(t) 


exp(-  *lt2) 


exp(-  Jit2 
exp(-  A  t2 


+  2/\igt1(t-t1)  ) 

+  2  Vg(t2+{l-g)tjt  -  2^1g(t2+t2-gt1t2) 


,  t  6  £o»t^) 

,  t  €  C^.tg) 
) ,  t  £  £t2 ,  t^ ) 


(2.4.12) 


1 


i 


§2.5  PROBABILITY  DENSITY  FUNCTIONS 


The  probability  density  function  f(t)  of  the  lifetime  T  is  related 
to  the  failure  rate  r(t)  and  the  reliability  function  R(t)  by  Eq. (2.2.7). 
Define 

f(t)  s  r(t)R(t)  >  regular  lifetime  p.d.f.  (2.5.1) 


and 


t6tVi'V 

f(t)  -  f  (t)  - 
n 

rn(t)R(t) 

(2.5.2) 

Thus, 

fx(t) 

3  r1(t)R1(t) 

.  *£[0,^) 

f2(t) 

3  r2(t)R1(t)R2(t) 

,  t^Ctlft2) 

f(t)  .  « 

• 

• 

• 

fn(t) 

• 

• 

• 

» 

n-1 

S  r„(t)Rn(«n  R j < > 

j»l 

*  t^ttn-l,tn) 

(2.5.3) 

Note  that 

tSCO.t^ 

+J 

0<M 

II 

-P 

H 

<H 

t 

(2.5.4) 

(1)  FIXED  REDUCTION  PROBABILITY  DENSITY  FUNCTION 

For  fixed  reduction  failure  rate,  Eq. (2.3.6),  (2.4.10)  and  (2.5.2) 

give 

f  (t)  -  (  r(t)  -  (n-1)^  )  R(t)  exp(  (n-1)  (t-J$nt. )  A  )  (2.5.5) 

n  x 

°  o 

»  (  f(t)  -  (n-1)^  R(t)  )  exp(  (n-l)(t-3^nt1)^  )  (2.5.51) 

«  -  2  &  M  2A*«tx 

fn(t)  »  £  2X*'t  -  (n-1)^  }  exp(  -  A*t2  +  (n-1)  (t-^rt^)^  ) 

and  the  shape  of  f(t)  is  sketched  in  Fig. 2. 5.1. 


(2.5.6> 


Fig.  2.5.1  Lifetime  probability  density  function  with  fixed  reduction 


(2)  PROPORTIONAL  REDUCTION  PROBABILITY  DENSITY  FUNCTION 


For  proportional  reduction  failure  rate,  we  have  for  0(  =  2  , 


’2  A*t.oxp(-  ^t2) 

2  .expj-  t2  +  2  gt^t-t^j 

f(t)«K  2  A  [t-g( t2+(  1-g) tx)] exp{- # [t2-2s( t2+(  1-g) t± ) t 

+2G(t2+t2rgt1t2)]  j> 

The  shape  of  f(t)  is  sketched  in  Fig.  2.5.2. 


,  t€  Co,tx) 

,  t£(tlft2) 


,  t^Ct2,t3) 


(2.5.7) 


Note  that  f^(t)  and  f2(t)  are  the  same  for  both  fixed  and  proportional 

reductions  for  a  given  failure  reduction  factor  g  (see  Fig.  2.5.1  and  2.5.2). 

Suppose  there  is  only  one  maintenance— cycle  and  a  large  operative— cycle  time. 

In  this  case  both  fixed  and  proportional  reductions  are  equivalent.  The 
o 

curves  of  f(t)  and  f(t)  (for  Ci  »  2)  are  illustrated  in  Fig.  2.5.3,  and 


f(t)  ■  2Alt  exp(— >ft2)  ,  t^O 

f  2  A*t  exp(- A*t2)  ,  t  £ [0,t  ) 


t  l (2)ft  -  /\)exp{-A*t  +  (t-tx)/^  ,  t^.^  ;  ^  3 

Let  t  be  the  point  that  both  curves  intersect,  that  is, 
f(t*)  -  f(tp)  *  0 


(2.5.8) 


**  (2.5.9) 


(2.5.10) 


[-[t'-tjA) 


Expanding  the  exponential,  we  have  approximately, 
t#  * 


*  tX  ~ 


for  large  t. 


(2.5.11) 


(2.5.12) 


°<t) 


(2.5.12') 


0 


I  ■  ■  * 

2/A, 


Fig.  2.5.2  Lifetime  probability  density  function  with  proportional  reduction. 


JU 


And 


*1  - 

00 

tw  i  tx 

(2.5.13) 

f(t)  > 

f(t) 

for 

t<t* 

(2.5.14) 

ct 

V 

0 

f(t) 

for 

t  >  te 

(2.5.15) 

Ue  conclude  that  for  large  operative-cycle  time  and  a  single  maintenance 
operation,  the  lifetime  probability  density  function  with  failure  reduction 
is  always  greater  than  the  regular  lifetime  p.d.f.  whenever  t  >t^  ,  namely. 


?  tCs 


f(t)  >  f(t) 


(2.5.16) 


Fig.  2.5.3  Asymptotic  behavior  of  lifetime  p.d.f.  for  a 
single  maintenance. 


2.6  OPERATIVE  TIME  AND  MAINTENANCE  TIME  DISTRIBUTIONS 


We  define  an  integer  random  variable  N^  to  be  the  number  of  Q-M  cycles 
before  failure.  It  is  obvious  from  Fig.  2.1.2  that 


n 


t  €  ft  ,t  ) 
n— l  n 


(2.6.1) 


for  n  »  1,2,3, .. .  ,  and  t^aO. 

Note  that  for  a  renewal  process  (Parzen  (1962))  the  expectation  of  is 

called  the  renewal  function.  We  shall  use  the  expected  value  of  N.  , 
namely,  E(  N^  }  ,  to  compute  the  approximate  mean  uptime  and  the  approximate 
mean  lifetime  for  a  periodic  operative-maintenance  schedule  (  Eq. (2.6.22) 
and  (2.6.23)  ). 


Since  the  system  may  fail  during  the 
th 


th 


O-cycle  (i.e.  in 


or  during  the  n  M-cycle  (i.e.  in  (W  >•  *•  further  define  the 


‘Vi'V  > 


following  two  randowm  variables 
lifetime  T  ; 


u)n 


and  T 


din 


related  to  the  system 


For  n  »  1,2,3, . . . 

T  ,  5  operative  time  before  failure  in  the  n^  O-cycle 

u,n  th 

(the  n  O-cycle  operative  time) 

T , ,  *  maintenance  time  before  failure  in  the  n^  M-cycle 

“In  . , 

(the  n  M-cycle  maintenance  time) 


such  that 


T  =  i 


t  +  T  , 
n-1  ujn  ,  if 


%i  +  Tdln  ,  if 


©-© 

©-© 


when  N  =n 


when  Nt=n 


(2.6.3) 


It  is  evident  that 
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0  <  T  ,  <  u 

v  uln  v  n 


(2.6.4) 


and 


Furthermore, 


0  s< 

T  <  d 

d|n  *  n 

(2.6.5) 

T  , 
u|n 

<un 

T,,  -  o 

d|n 

(2.6.6) 

T 

din 

S<dn 

T  .  m  U 
u|n  n 

(2.6.7) 

The  relationship  of  T  ,  Tujn  and  Td|n  i^us'fcra'ted  ^  FiS*  2.6.1. 


System 

lifetime 


t 


H - b 


-{h 


u|n 


d|n 


or 


n-1 


n 


Possible  failure  instant 


Fig.  2.6.1  Lifetimes  relationship  before  failure 


and 


Next,  we  shall  compute  the  probability  density  functions  of 
,  denoted  respectively  by 


f™  and 

u|n 


T 


u|n 


Suppose  =  n  and  the  p.d.f,  of  T  assuues  the  shape  shown  in 
Fig.  2.6.2  in  which  the  ?n(t)  is  defined  by  Eq.(2.5.2)-(2.5.3). 


f  (t) 
n 


rT  (t) 

d|n 


Fi«.  2.6.3  P.d.f.  Of  TuJn 


Fig.  2.6.4  P.d.f.  of  . 


•H.J 


(t) 


u|n 


.  (t) 

d|n 


1  v  t  , 
n  n-1 

U  +  D 
n  n 


U  +  D 
n  n 


V 


u 

n 

U  +  D 
n  n 

f  (*  +t> 

n  n _ 

U  +  D 


n  n 

0 


,  0  $  t  <  un 

,  t  =  u 

n 

,  otherwise 

.  t  -  0 

•  0  <  *  d„ 

,  otherwise 


(2.C.13) 


(2.6.14) 


Furthermore,  from  Eq.(2.5.3), 


n 


D 

n 


n 


fn(x)  dx 


J  t, 


n-1 


n-1 


"  rn(x)Vx)  TT  VV 


n-1 

n-1 


Hvv 


j-i 


j-i 

[  1  -  w  ) 


•  • 

J  ' 


f  (x)  dx 
n 


(tn  n-1 

r  rn(x)Vx)  Tl  W  dx 


j=l 


ft  MV  (",WVl 


(2.6.15) 


! 


I 

(2.6.16)  ! 


I 


£4.  (2.6.10)  and  (2.6.16)  give 
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n  w 
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In  particular, 


n  *  1 


p{Nt  .  i}  .  1  - 


The  average  number  of  0-M  cycles  is  given  by 
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th 

Similarly,  the  k  moment  of  N  is  given  by 

E  (Ntk)  5  IZ  “k-pfNt-n! 
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n»l  U  j-l 
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(2.6.17) 


(2.6.18) 


] 


(2.6.19) 


(2.6.19*) 


(2.6.20) 
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Note  that  If  there  is  no  maintenance  operation,  or  the  operative-cycle 
time  becomes  infinite,  then 

t,  -*  oo  n  (t  )  -»  0 

1  11 

and  Eq. (2.6.19)  becomes 

lim  EC  Nt  )  »  1  (2.6.21) 

t^-*  QO 

That  is,  there  is  only  one  0-cycle.  (See  Table  4.4.1  and  Fig.  4.4.3  for 
a  numerical  illustration) . 

We  may  use  the  mean  0-M  cycles  E[  N..  3  to  compute  the  approximate 
mean  uptime  and  the  approximate  mean  lifetime  for  a  periodic  operative- 
maintenance  schedule  (  (2.1.1)  ). 

Define 

Approximate  mean  uptime  *  u. (  EC  \  )  (2.6.22) 

Approximate  mean  lifetime  *  (  u  +  d  ).(  eCNj.3  -  %  )  (2.6.23) 

We  shall  see  that  Eq. (2.6.22)  and  (2.6.23)  give  consistent  results  as 
— ao  (See  Table  4.4.1  and  Fig.  4.4.4-4.4.5) . 


u 


§2.7  REPAIR  TIME  DISTRIBUTION 


Whan  the  system  fails,  it  i3  in  the  stste  for  a  complete  repair, 

replacement  or  renewal.  The  time  spent  is  a  random  variable  called  the 
repair  time  Tp  with  distribution  function 

^P(t) 

and  probability  density  function 

p(t) 

For  simplicity,  Tp  is  assumed  to  statistically  independent  of  the  system 
lifetime  T  ,  and  Tp  is  exponentially  distributed. 


p(t) 

*  p  exp(  -  pt) 

,  t  JO 

(2.7.1) 

ct 

a  1  -  exp (  -  pt) 

,  tJO 

(2.7.2) 

We  shall  see  that  other  distributions  may  be  used  and  the  independence 
assumption  may  be  dropped,  but  we  have  to  handle  double  integrals,  instead 
of  single  integrals  in  evaluating  the  probability  of  stochastic  availability 
(see  §  3.2) . 


CHAPTER  3 


STOCHASTIC  AVAILABILITY 


§3.1  DEFINITIONS  OF  STOCHASTIC  AVAILABILITIES 

We  define  the  following  random  variables  with  reference  to  Fig.  2.1.2 
and  2.6.1  : 

*  total  uptime  before  failure  when  the  system  is  operative 
5  total  downtime  during  which  the  system  is  not  operative 
a  \  u.  +  T  , 

J“1  n-  1,2,3,... 

n-1 

*  dj  +  Td|n  +  TF 

where  the  summations  are  defined  to  be  zero  for  n  ■  1. 

The  following  stochastic  availabilities  (A)  are  defined  : 

(1)  Stochastic  Cycle  Availability  (Stochastic  Availability  for  a  complete 
Renewal  Cycle) 


T 

up 

T 

down 

T  , 
up  |n 

T 

down|n 


T  +  T 
up  down 


(3.1.1) 
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(2)  Stochastic  Availability  for  the  n  0-M  Cycle  A„  (  A  ) 

■"  ■  N^.*n  n 


.  .*  ,  - seifi - 

V"  n  t  *  T 

up | n  down | n 


n-1 

Eu.  +  T  , 

J  u|n 
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>  u.  +  T  ,  +  \  d.  +  T. ,  +  T_ 

j  ujn  / _  j  d|n  F 


Et 

j-i 


(3.1.2) 


(3.1.2*) 


(3)  Stochastic  Availability  for  the  nm  O-cycle  A. 


n-1 

Eu.  +  I  . 
j  u|n 

J-1 


0  n-1 


n  Euj*vEvtf 


(3.1.3) 


T  mO 
d|n 


(3.1.3*) 


(4)  Stochastic  Availability  for  the  n  14-cycle  A^j 


A,.  *  — 

A  n 
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(3.1.4) 


(3.1.4*) 
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( 5)  Stochastic  Availability  at  a  Finite  Time  in  the  11th  0-M  Cycle  t 


(Finite  Time  Stochastic  Availability) 
r  n-l 


.t  , 


21  uj  *  "in<t'Tui») 


“j  ♦Tr 


n-l  n-l 
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l  J-1  J-1 


t-r  +t  <  t 

n  n 


(3.1.5) 


Hence, 


.t-r 

An  n 


.t-t 
An  n 


(3.1.5') 

(3.1.5") 


(6)  Stochastic  Availability  for  an  Age  Replacement  Time  A 


AR 


t  uptime  before  failure  in  (Q,t_) 

A  R  a  - - - JL_ 

AR 


uptime  before  failure 

in  C0,tR) 


+  Total  downtime 


(3.1.6) 


where  tR  is  the  age  replacement  time.  See  §4.2  for  more  detail. 


We  note  that  the  stochastic  availabilities  defined  above  are  in  the  form 

uptime 

A  "  uptime  +  downtime  (3.1.7) 

where  the  uptime  means  that  the  system  is  operative  and  the  downtime  means 
that  it  is  not  operative  such  as  in  the  maintenance  or  failure  state. 

Since  the  uptime  and/or  the  downtime  are  in  general,  nonnegative  random 
variables,  the  stochastic  availability  is  therefore  random  and  less  than  unity. 


§3.2  PROBABILITY  DISTRIBUTIONS  OF  STOCHASTIC  AVAI LI ABILITY 


We  have  defined  various  stochastic  availabilities  A  in  the  provious 
section.  We  shall  derive  the  probability  of  the  form 

p{A^l-t}  ,  0<|<1 

in  this  section  and  in  §  4.2.  We  shall  begin  with  the  stochastic  cycle 
availability. 

By  the  total  probability  theorem, 
p  {  Ac*1  -  *} 

•  fz  pf  AC*1'£  I  v } .  p(vn} 

nail 
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(3.2.2) 


Consider 


A.  Jl-! 


Eu.  +  T  . 

j  ujn 

3*1 _ 

n-1  n-1 

Eu.  +  T  ,  +  N  d.  +  T_ 

j  u|n  j  F 


2  1  -  f 
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(3.2.3) 


A  J 


The  domain  of  definition  of  T^n  and  Tp  is  shown  in  Fig.  3.2.1. 

Let  f  «  (...)  be  the  joint  probability  density  function  of  T  .  and  T_ 
Tu|n’TF  u‘n  F 

By  the  independence  assumption  (§  2.7),  we  have 
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Let  z  =  x  +  t^_^  and  from  Eq. (2.0.13),  we  have 

p£->  n_1 
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Next,  consider 
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The  domain  of  definition  of  T, .  and  T_  is  shown  in  Fig.  3.2.2. 

d  |  n  r 

Let  f„  _(.,.)  be  the  joint  probability  density  function  of 
d|n*  F 

Tdjn  and  Tp  and  by  the  independence  assumption  (§  2.7), 
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Fig.  3.2.2  Domain  of  definition  of  T.. 
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Then, 
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Let  z  *  x  +  Vn  and  from  Eq.(2.6.14),  we  have 
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From  Eq. (3.2.1),  (3.2.2),  (3.2.10),  (3.2.16),  (2.6.11)  and  (2.6.12),  we 


have 
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Eq.(3.2.17)  and  (2.6.10)  give  the  probability  distribution  of  the 
stochastic  cycle  availability  : 
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The  finite  tine  stochastic  availability  can  be  computed  similarly  by 


the  following  changes  (see  Eq. (3.1.5*)  and  (3.1.5"))  : 
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Then, 
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The  probability  distribution  pf  the  stochastic  availability  for  an 
age  replacement  time  will  be  derived  in  §  4.2  (Eq. (4.2.9) ) . 


The  probability  distributions  are  summarized  in  the  following  : 


PROBABILITY  DISTRIBUTIONS 
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g  3.3  ASYMPTOTIC  DISTRIBUTION 


SO 


We  shall  study  the  behavior  of  p[acJ1-e|  as  the  operative-cycle  time 
increases  indefinitely  (that  is,  no  maintenance,  — ►  ®  ) .  This 

asymptotic  value  will  be  denoted  by 
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Recall  from  Eq.(3.2.18), 
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is  the  Laplace  transform,  or 
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Some  asymptotic  values  are  tabulated  in  Table  3.3.1. 
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0.1 
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0.23878 

0.7063 

0.15 

0.1 

1 

15/17 

0.21227 

0.7231 

Table  3.3.1  Some  asymptotic  probability  values  of 
stochastic  cycle  availability. 


(3.3.13) 


CHAPTER  4 


OPTIMUM  SYSTEM  DESIGN 


§4.1  PERIODIC  OPERATIVE-MAINTENANCE  POLICY 

We  have  mentioned  in  §  2.3  that  for  the  case  of  system  performance  and 
management,  it  is  more  efficient  to  have  a  periodic  operative-maintenance 
schedule,  namely,  constant  operative-cycle  time  u  and  constant  maintenance- 
cycle  time  d  (Eq. (2.1.1) ) .  Normally,  the  maintenance-cycle  time  d  is 
more  or  less  fixed  by  the  system  parameters  and  we  may  free  to  choose  an 
optimum  operative-cycle  time  u*  such  that  certain  performance  index  is 
satisfied  in  order  to  give  an  optimum  system. 

Classically,  as  a  performance  index,  the  minimization  of  a  long  term 
average  operating  cost  or  the  maximization  of  the  system  availability  has 
been  proposed  (Barlow  and  Proschan  (1965)).  Since  we  have  defined 
availability  as  a  random  varoable,  we  may  choose  the  criterion  of  maximizing 
the  probability  of  stochastic  availability  with  respect  to  the  operative- 
cycle  time  to  be  an  optimum  policy.  In  this  case  we  shall  have  a 
probabilistic  guarantee  instead  of  using  averaged  quantities.  We  shall,  see 
in  Chapter  5  that  this  concept  is  equivalent  to  that  of  minimizing  an 
appropriate  system  cost  function. 

If  we  plot  p{ac^1  -£ j1  versus  u  where  Ac  is  the  stochastic  cycle 
availability  and  1,  we  expect  the  following  : 

(1)  If  there  is  no  failure  reduction,  P|AcJl-t}  increases  with  u  because 
the  more  frequent  the  number  of  preventive  maintenance,  the  less  is  the 
availability. 

(2)  If  there  is  a  failure  reduction  due  to  the  maintenance  operations,  there 
may  exist  an  u  such  that  P  |  A_  £  1— £ }  is  a  maximum. 

(3)  As  u  increases  indefinitely,  there  is  essentially  one  operative  cycle 
and  no  maintenance  operation,  therefore,  (see  §  3.3) 

lim  p{Ac}l-l(  -  Pffl 
u-»co 
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The  above  reasoning  is  illustrated  in  Fig.  4.1.1  and  the  conjecture 
will  justified  by  examples  in  |  4.4  • 


?{ V'  *-‘l 


Fig.  4.1.1  Probability  of  stochastic  cycle  availability  versus 

constant  operative-cycle  time  for  periodic  0-M  schedule. 


From  Fig.  4.1.1,  v/e  may  choose  an  optimum  operative-cycle  time  u  in 
one  of  the  following  criteria  : 


(a)  If  there  a  failure  reduction,  choose 
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is  achieved. 
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(b) 
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is  a  failure  reduction. 
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for  u^u 

where 
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is  no  failure  reduction, 
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r  (4.1.2) 


such  that  the  total  operation  cost  is  as  small  as  possible,  namely, 
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00 


>  (4.1.3) 


The  above  criteria  apply  the  concept  of  maximizing  the  probability  of 
stochastic  cycle  availability  and  probabilistic  inequalities  in  system 
design.  Other  stochastic  availabilities  may  be  used.  In  particular, 
the  stochastic  availability  for  an  age  replacement  time  will  be  considered 
in  §  4.2  and  the  concepts  will  illustrated  numerically  in  §  4.4.  We 
shall  show  that  these  concepts  are  equivalent  to  those  of  minimizing  the 
probability  of  an  appropriate  cost  function  and  the  corresponding 
probabilistic  inequalities  involving  this  cost  function  in  Chapter  5. 


3  4.2  AGE  REPLACEMENT  POLICY 


00 


When  the  sytem  failure  rate  increases  with  age  and  system  failure  is 
costly,  it  is  usually  scheduled  to  replace  or  renewal  the  entire  system 
before  it  has  aged  too  greatly.  This  is  called  an  ane  replacement  policy 
(Carlow  and  Proschan  (1905)}  which  replaces  the  system  at  a  time  after 

its  installation  or  at  failure,  whichever  occurs  first.  The  time  t_  is 

n 

called  the  age  replacement  time.  We  shall  find  an  age  replacement  time  for 
the  repairable  system  with  an  age  and  maintenance  dependent  failure  rate 
using  the  concept  of  stochastic  availability. 


Suppose  the  age  replacement  time  t_  is  chosen  to  be  at  the  end  of  the 

K  operative-cycle  (Fig.  4.2.1),  that  is,  the  system  is  replaced  instead  of 

til 

doing  preventive  maintenance  after  the  a  1  0-cycle,  where  the  value  of  K  , 
or  equivalently,  tR  is  to  be  determined.  We  have,  from  Fig.  4.2.1, 
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Fig.  4.2.1  An  age  replacement  schedule 
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(4.2. 7') 


where 


^  (t)  «  1  -  exp  (-  jit)  ,  t  ^  0 

is  the  distribution  function  of  the  repair  time  (Eq. (2.7.2) ). 


Note  that 


D  =  0 
K 


(4.2.0) 


Now  we  compute  the  probability  of  the  stochastic  availability  for  an  age 
replacement  time.  Oy  the  total  probability  theorem,  Eq. (4.2.5)-(4.2.G) , 
and  the  probability  distributions  summary  in  §  3.2, 

p  { >  1  - £  1 

■  E  p{Al!Ul-£!  "t-- W V"}  *  flp{A^1-£  I  v'}pK™| 
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E  .  E  S^n - n v« u» * v 

n»l  U  +  D  n-K+1  6j»l  J  j-l  J 
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“I  am  >'  l-£l 

K— 1  ,,  A,  n  m  _  K-l  K  K-l 

r’4^-o,  •  -gw, 
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As  t  t  cd  ,  i.e.  K  -*  a>  ,  we  have 

AAR  “  Ac 

V  ® 


(4.2.9) 


(4.2.10) 


and 


*fAjS>l*Sl  -  pfAo>1-‘l 


(4.2.11) 


For  a  periodic  operative— maintenance  schedule,  we  expect  that  the  graph 

of  P  j  A  n  X 1-  £  [  versus  t„  to  have  the  3hape3  shown  in  Fig.  4.2.2  for 
[  AR  I  R 

system  with  failure  reduction.  The  curves  are  justified  in  §  4.4.  These 
curves  are  similar  to  those  in  Fig.  4.1.1.  Hence  we  propose  the  following 
age  replacement  policy  : 

(A)  If  the  plot  of  P  {  }l-£  | 

maximum,  then  choose 


versus  t  for  a  given  u  shows  a 

K 


where 


t  a  t 
R  R 


max 
t. 

is  achieved. 

(0)  If  the  plot  of  P 


*  ?  {  aar  ^  i  1 


Y  (4.2.12) 


versus  t  for  a  given  u  shows  a 
K 


maximum,  then  choose  t  such  that 

•  A 


where 


p{AI? 

l  1  -tj 

^  1 

- s?  for  s  *  *i 

*  * 

1  -  ( pe 

,  8>o 

p<  5 

-} 

(for  the  given  u  ) 
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f  (4.2.13) 
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(C)  if  the  plot  versus  tR  shows  no  peak,  thon  chooso 

t_  to  be  the  first  time  such  that 


P  (  ^  1  y  1  -  S 

I  Alt  '  1 1  '  Q 

v/h ere 

fQ  s  1  -  (  Pe  - J)  ,  J>o 


(4.2 


In  any  of  the  above  cases, 
K^1  O-cycle. 


t_  is  chosen  to  be  at  the  end  of  a 


Fig.  4.2.2  Probability  of  age  replacement  stochastic  availability 
versus  age  replacement  time  for  failure  reduction  case 


ul 

Other  types  of  ago  replacement  policies  can  be  proposed.  Usually, 
system  failure  during  operation  is  a  serious  problem  than  when  it  fails 
during  maintenance.  A  simple  age  replacement  policy  based  on  this  idea 
is  stated  below  ; 

The  repairable  system  with  maintenance  schedule  is  replaced 
when  the  probability  of  system  failure  during  operation  is  q  , 
where  q£(0,l),  say,  q  =  Vi  . 


We  have  from  Fig.  2.G.1,  the  probability  distributions  summary  in  §  3.2, 
and  the  total  probability  theorem, 

P  {  Failure  during  operation  | 

a  y  P  |  Failure  in  the  nth  0-M  cycle  J  Nt*n  |  P  {  N^=n } 
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(4.2.15') 


For  an  infinite  operative-cycle  time  (no  maintenance),  u  oo, 
oo  ,as 


lia  II  U„  .  Jn 

u-»  oo  n»l  *  0 


f(z)dz 


(4.2. 16) 


That  is,  the  system  will  ultimately  fail  as  expected.  With  the  age 
replacement  time  t^  chosen  to  be  the  end  of  an  operative-cycle,  the 
probability  curve  P  f Failure  during  operation^  versus  tR  is  a 
monotonic  increasing  function  (Eq.  (4.2.15-) ) ,  and  is  depicted  in  Fig.  4.2.3. 


§  4.3  COMPUTATION 


We  have  studied  a  repairable  system  with  age  and  maintenance  failure  rates, 
stochastic  availabilities,  periodic  operative-maintenance  policy,  and  age 
replacement  policy.  We  shall  study  various  characteristics  of  the  repairable 
system  via  exasqples. 

In  order  to  evaluate  various  probabilities  of  stochastic  availabilities, 
a  computer  program  is  developed  for  computational  purpose.  The  program  is 
written  in  PASCAL  language  and  is  simulated  by  an  APPLE  II  PLUS  (1080  model) 
microcomputer  with  48K  memory. 

normally  we  assume  that  the  statistics  of  the  system  failure  rate 
function  are  known  or  can  be  estimated.  We  shall  mainly  concern  with  the 
periodic  operative-maintenance  policy,  that  is,  fixed  maintenance-cycle  time 
(d)  and  constant  operative-cycle  time  (u)  .  Variable  operative-  or 
maintenance-  cycle  times  can  be  considered  by  a  slight  modification  of  the 
computer  program. 

The  program  has  to  evaluate  the  probability  distributions  of  the 
stochastic  availabilities.  In  particular,  the  stochastic  cycle  availability 
(Eq.(3.2.13))  involves  an  infinite  summation.  An  upper  limit  of  the 
summation  is  assigned  to  be  an  input  to  the  program  and  the  program  checks 
the  difference  between  two  consecutive  partial  sums  whether  it  is  less  than 
10  which  is  the  accuracy  of  the  machine.  Simpson's  rule  is  used  to 
evaluate  all  the  integrals  with  the  number  of  integration  points  as  an  input. 

The  number  of  integration  points  is  about  30-50 'in  the  examples  (5  4.4). 

The  number  of  summation  varies  from  10  to  1000  or  more  and  it  depends  on  the 
failure  reduction  type  and  failure  reduction  factor  3  .  Since  all 
quantities  involved  in  computation  are  nonnegative,  various  checks  are  made 
to  ensure  this  nonnegativity.  Underflow  or  division  by  zero  i3  suppressed 
to  give  correct  results.  Since  the  machine  has  limited  accuracy  and  low 
speed,  numerical  difficulties  occur  quite  frequently,  if  the  program  were 
Witten  for  a  large  computer  with  double  precision,  accuracy  and  speed  will 
be  improved  but  with  a  tradeoff  of  higher  cost. 


I 


The  inputs  to  the  computer  program  are  : 


(1)  £  ,  0<  t<  l 

(2)  0(  ,  usually  or  a  2 

(3)  ^1  ,  Agj  »  where  Afl  -  -  A^expMn-D/lO) ,  n»l,2 

■  if  A.  *  A  ,  then  A  ■  A  =  constant 
1  qd  n 

(4)  Reduction  type  : 

0  >  no  reduction 

1  >  fixed  reduction 

2  a  proportional  reduction 

( 5)  Reduction  factor  :  0  4  S 4 1 

(6)  Maintenance-cycle  time  :  d 

(7)  Operative-cycle  time  :  u 
(3)  Maximum  number  of  summation 
(9)  Integration  points,  e.g.  30. 

The  outputs  from  the  program  are  : 

(1)  p{Ac*l-t} 

(2)  EL*  N.  1 

(3)  t 

(4)  p{a^^1-£ 

(5)  P  |Aq  £  l  -t 

(6)  P  {Ajj  *  l  -«} 

(7)  p|a^1-(| 

(8)  P  {  -  n  [ 

(9)  Approximate  mean  uptime 

(10)  Approximate  mean  lifetime 

A  numerical  illustration  of  the  computer  program  is  given  below. 

A  complete  listing  of  the  program  is  given  in  Appendix  2. 


DATE  :  APRIl  Hi  1932 


**  STOCHASTIC  AVAILABILITY  ** 
<**  FIXED  REDUCTION  **> 


REDUCTION 

FACTOR  = 

0.50000 

c  - 

0.10 

APLHA  = 

O 

4m 

LAMBA  = 

0.1000 

MU  = 

1.0000 

MAINTENANCE-CYCLE  TIME  = 

0.01000 

INTEGRATION  =  36  POINTS 


OPERATIVE-CYCLE  TINE  =  2.00000 


N 

PCAOl-El 

EC  N(  T )  3 

ART 

PC AAR 3 

PC  AO  3 

PC  AM] 

PC  AN] 

PCN(  T)*M] 

i 

* 

0.005281 

0.04 

2.000 

0.196705 

0.134662 

0.0013994 

0.1346804 

3.95953E-2 

4» 

0.026384 

0.19 

4.010 

0.338639 

0.235096 

0.0024982 

0.2851142 

7.45501E-2 

3 

0.068266 

0.49 

6.020 

0.442228 

0.416777 

0.0029526 

0.4167966 

1.0U17E-1 

4 

0.129460 

0.96 

9.030 

0.514615 

0.525566 

0.0032838 

0.5255858 

1.17100E-1 

5 

0.204083 

1.57 

10.040 

0.563100 

0.614410 

0.0035316 

0.6144296 

1.22117E-1 

6 

0.284313 

2.28 

12.050 

0.594136 

0.686677 

0.0037008 

0.6866963 

1.17430E-1 

7 

0.362529 

3.01 

14.060 

0.613240 

0.745359 

0.0038123 

0.7453783 

1.05454E-1 

3- 

0.432853 

3.73 

16.070 

0.624396 

0.792967 

0.0038789 

0.7929869 

8.91043E-2 

9 

0.491732 

4.37 

18.080 

0.630594 

0.831574 

0.0039111 

0.8315921 

7.11892E-2 

10 

0.538133 

4.91 

20.090 

0.633874 

0.862372 

0.0039147 

0.8628907 

5.39567E-2 

11 

0.572528 

5.33 

22.100 

0.635509 

0.888245 

0.0038982 

0.8882622 

3.89887E-2 

12 

0.596692 

5.66 

24.110 

0.636271 

0.908812 

0.0038655 

0.9083285 

2.67000E-2 

13 

0.612809 

5.88 

26.120 

0.636599 

0.925485 

0.0033200 

0.9255010 

1.74357E-2 

1* 

0.623035 

6.04 

23.130 

0.636724 

0.939002 

0.0037667 

0.9390177 

1.09343E-2 

15 

0.629219 

6.13 

30.140 

0.636764 

0.949965 

0.0037052 

0.9499302 

6.53417E-3 

16 

0.632786 

6.19 

32.150 

0.636772 

0.958857 

0.0036409 

0.9583719 

3.73393E-3 

17 

0.634751 

6.23 

34.160 

0.636770 

0.966072 

0.0035703 

0.9660353 

2.04153E-3 

13 

0.635785 

6.25 

36.170 

0.636767 

0.971926 

0.0034987 

0.9719337 

1.06844E-3 

19 

0.636307 

6.26 

38.180 

0.636764 

0.976683 

0.0034257 

0.9766957 

5.35445E-4 

20 

0.636558 

6.26 

40.190 

0.636762 

0.980547 

0.0033518 

0.9805588 

2.57031E-4 

21 

0.636674 

6.26 

42.200 

0.636762 

0.983689 

0.0032777 

0.9837007 

1 . 18218E-4 

22 

0.636725 

6.27 

44.210 

0.636761 

0.986246 

0.0032033 

0.9862578 

5.21078E-5 

23 

0.636747 

6.27 

46.220 

0.636761 

0.988329 

0.0031303 

0.9883393 

2.20159E-5 

24 

0.636755 

6.27 

48.230 

0.636761 

0.990033 

0.0030537 

0.9900427 

3.91778E-6 

25 

0.636759 

6.27 

50.240 

0.636761 

0.991418 

0.0029856 

0.9914276 

3.46367E-6 

26 

0.636760 

6.27 

52.250 

0.636761 

0.992559 

0.0029136 

0.9925676 

1 .2901 IE-6 

27 

0.636760 

6.27 

54.260 

0.636761 

0.993488 

0.0028449 

0.9934969 

4.60982E-7 

29 

0.636761 

6.27 

56.270 

0.636761 

0.994259 

0.0027774 

0.9942664 

1.57929E-7 

29 

0.636761 

6.27 

58.280 

0.636761 

0.994889 

0.0027089 

0.9943979 

5.19149E-8 

OPERATIVE-CYCLE  TINE  =  2.00000 
PC  SCA>  0.90  3  =  0.636761*  C29] 
NEAN  0-M  CYCLES  =  6.26696 
MEAN  UPTINE  =  1.15339E1 
NEAN  LIFETIME  =  1.15916E1 


oo 


Table  4.3.1 


Numerical  illustration  of  the  PASCAL  computer  program 


4.4  NUMERICAL  EXAMPLES 


In  this  section  examples  will  be  given  to  illustrate  various  concepts 
developed  so  far  and  to  explore  the  characteristics  of  the  repairable  system 
with  age  and  maintenance  dependent  failure  rates.  In  particular,  Fig.  4.1.1 
and  4.2.2  are  verified.  The  asymptotic  values  in  Table  3.3.1  will  be  used 
in  all  the  examples. 


EXAMPLE  4.4.1 


We  consider  the  repairable  system  with  various  failure  reduction 
criteria,  namely, 

(a)  no  reduction  (NR) 

(b)  fixed  reduction  (FR) 

(c)  proportional  reduction  (PR) 

for  a  piecewisely  linear  failure  rate  with  the  following  parameters  : 

<X  =»  2,  0.1,  ^  *  0.1,  /i  »  1,  d  =«  0.01,  g  3  0.5 

The  failure  rate  functions  for  fixed  and  proportional  reduction  are  plotted 
in  Fig.  4.4.1.  The  values  of  P  {  A  £  0.9}  ,  EC  N.  J,  approximate  mean 
uptime  (Eq. (2.6.22) )  and  approximate  mean  lifetime  (Eq. (2.6.23) )  for 
(a)  NR,  (b)  FR,  and  (c)  PR  are  computed  using  the  PASCAL  computer  program 
(Appendix  2)  for  various  operative-cycle  time  u  ,  and  they  are  tabulated 
in  Table  4.4.1. 

P  ■[ Ac  0 . 9  }•  is  an  increasing  function  of  u  for  no  reduction  and  it 
has  a  maximum  for  both  types  of  failure  reduction.  For  all  types  of  failure 
rates  P  f  A  £  0.9  f  approach  the  asymptotic  value  P  (=»  0.5794,  see 

C  *  CD 

Table  3.3.1).  These  curves  are  plotted  in  Fig.  4.4.2  which  when  compared 
with  Fig.  4.1.1,  the  statements  made. in  §  4.1.1  are  verified. 

From  appendix  1,  Eq.(A1.6),  the  regular  no  failure  reduction  Weibull 


mean  lifetime  is 
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4.4.1  Linear  failure  rate  for  fixed  and  proportioned  reduction 


Approximate  mean  lifetime  Approximate  mean  uptime 


EXAijPLE  4.4.2 


-e  shall  study  the  characteristic  of  parameter  variations  of  the 
repairable  system  with  fixed  reduction  (proportional  reduction  is  similar) . 
The  values  of  p{  l-i|  versus  the  operative-cycle  time  u  for  of  =  2 
(fixed  reduction)  are  computed  by  the  PASCAL  computer  program  and  they  arc 
tabulated  in  Table  4.4.2  for  different  variations  of  parameters  : 

i  ^  ,  p  ,  d  ,  g 

The  asymptotic  values  P  of  the  stochastic  cycle  availabilities  have 
been  computed  in  Table  3.3.1.  The  corresponding  curves  are  plotted. 

(1)  REDUCTION  FACTOR  VARIATION  (s) 

Fig.  4.4.6  depicts  the  variations  of  reduction  factors  g  =  0,  0.5, 

and  0.9  (columns  (1),  (2),  and  (3)  in  Table  4.4.2).  The  larger  the 

reduction  factor  the  more  frequency  of  preventive  maintenance  is  favored 

* 

to  improve  the  system.  Thus  the  optimum  operative-cycle  time  u  is 
shorter  for  larger  g  to  have  more  maintenance-cycles.  When  there  is 
no  failure  reduction,  that  is,  g  *  0,  u*  approaches  infinite. 

Mote  that  all  curves  approach  asymptotically  P  ^  . 

(2)  :WVi;?TEr:AHCE-CYCLE  TIME  VARIATION  (d) 

The  variations  of  maintenance-cycle  time  (d)  is  illustrated  in 

Fig.  4.4.7  (columns  (2)  and  (4)  in  Table  4.4.2).  The  longer  the  u  , 

the  less  availability  of  the  system,  hence  it  favors  for  a  longer 

,  * 

optimum  operative-cycle  time  u 


p{ac^  1  -<i 

(o<=  2)  Fixed  reduction 

:  0.1 

■SBB 
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0.1 

0.15 

:  0.1 
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0.1 

0.1 

0.05 

.osf.io 

0.1 

0.1 

0.1 

u 

u:  1 

1 

1 

1 

1 

1 

1.5 

2 

1 

d:  0.01 

0.01 

0.05 

El 

0.01 

0.01 

0.01 

mm 

g:  0 

in 

e 

o 

Btfl 

0.5 

m 

0.5 

0.5 

0.5 

ESI 

• 

• 

• 

0.5 

0.4930 

0.5899 

:  0.7787 

0. 1139 

0.7827 

0.6256 

0.7078 

0.7797 

0.7424 

1 

0.5379 

0.6287 

0.7921 

0.4386 

0.8195 

0.7152 

0.7445 

0.3132 

0.7675 

1.5 

0.5527 

0.6363 

0.7813 

0.5189 

0.8289 

0.7519 

0.7516 

0.8197 

0.7714 

2 

0.5600 

0.6368 

0.7663 

0.5528 

0.8313 

0.7719 

0.7520 

0.8201 

0.7702 

2.5 

0.5644 

0.6346 

0.7505 

0.5701 

0.8324 

0.7841 

0.7500 

0.8185 

0.7674 

3 

0.5673 

0.6314 

0.7350 

0.5796 

0.8318 

0.7920 

0.7470 

0.8160 

0.7639 

3.5 

0.5694 

0.6277 

0.7203 

0.5849 

0.8305 

0.7973 

0.7437 

0.8131 

0.7602 

4 

0.5709 

0.6233 

0.7065 

0.5878 

0.8290 

0.8009 

0.7402 

0.0103 

0.7565 

4.5 

0.5721 

0.6199 

0.C93C 

0.5892 

0.8272 

0.6032 

0.7360 

0.0074 

0.7529 

5 

0.5731 

0.6162 

0.6017 

0.5896 

0.8254 

0.8047 

0.7334 

0.3048 

0.7495 

5.5 

0.5739 

0.6126 

0.6708 

0.5894 

0.8235 

0.8056 

0.7303 

0.8023 

0.7464 

6 

0.5745 

0.6091 

0.6607 

0.5888 

0.8216 

0.8060 

0.7274 

0.3000 

0.7434 

G.5 

0.5751 

0.6059 

0.6515 

0.5880 

0.8197 

0.8061 

0.7247 

0.7979 

0.7407 

7 

0.5756 

0.6029 

0.6431 

0.5871 

0.8178 

0.8059 

0.7222 

0.7960 

0.7383 

7.5 

0.5760 

0.6002 

0.6355 

0.5861 

0.8160 

0.8055 

0.7199 

0.7943 

0.7360 

3 

0.5763 

0.5976 

0.6286 

0.5851 

0.8142 

0.8050 

0.7179 

0.7929 

0.7340 

3.5 

0.5766 

0.5953 

0.6224 

0.5842 

0.3125 

0.8043 

0.7161 

0.7916 

0.7323 

9 

0.5769 

0.5932 

0.6163 

0.5833 

0.8108 

0.8036 

0.7145 

0.7905 

0.7307 

9.5 

0.5772 

0.5913 

0.6118 

0.5824 

0.8092 

0.8029 

0.7131 

0.7096 

0.7293 

10 

0.5774 

0.5896 

0.6073 

0.5817 

0.8077 

0.8021 

0.7118 

0.7888 

0.7281 

• 

• 

• 

i 

* 

* 

* 

1 

* 

; 

1 

V 

CD 

0.5794 

0.5794 

0.5794 

0.5794 

0.7859 

0.7859 

0.7063 

0.7859 

0.7231 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(0) 

(9) 

Table  4.4.2  P^A^l-f}  versus  u  with  variations  of  parameters. 


0(  a  2 


£  *  0.1 
A  *  0.1 


M  m  1 

g  *  0.5  (fixed  reduction) 


■« - - - - i - 1 . .  . .  . . .  .  -  u 

3  456  789  10 


versus  u  for  various  maintenance-cycle  time: 


(3)  VARIATION  OF  FAILURE  RATE  PARAMETER  (A) 

The  failure  rate  paramter  A  is  assumed  to  change  after  each  preventive 
maintenance  and  it  takes  the  form  (other  forms  may  be  proposed)  : 

A  a  A  —  (  ^  —  A.)  ®xp(~  i n  y  *  '  n»i,2f . . .  (4.4.4) 

n  oo  oo  l  xo 

The  failure  rate  function  r(t)  with  this  maintenance-cycle  dependent 
failure  rate  paramter  A  is  plotted  in  Fig.  4.4.0  (for  fixed  reduction). 
When  A  varies  from  A,  to  X  ,  we  shall  write 

0  X  QD 

^1  t  ^OD 

(as  0.05  f  0.1  in  Table  4.4.2  Column  (6)).  The  curves  of  P  | 

versus  u  for  X  »  0.05,  0.1  and  0.05  t  0.1  are  plotted  in  Fig.  4.4.9 

(columns  (5),  (2),  and  (6)  in  Table  4.4.2).  The  larger  the  A*  the 

higher  is  the  failure  rate,  hence  it  favors  for  more  failure  reduction,  or 

a  shorter  optimum  operative-cycle  time  to  have  more  14-cycles.  Thus 

u  is  shorter  for  A  *  0.1  than  that  for  A  ■  0.05.  For  the  maintenance 

dependent  A  >  as  u  is  small,  more  0-M  cyles  will  occur,  and  hence  the 

system  is  closer  to  that  for  X  »  A  .  As  u  ao  ,  less  number  of 

a? 

0-M  cycles  will  take  place  and  the  system  resembles  the  one  with  A-  ^  . 
Asynqptotically ,  the  values  P  =  0.5794  and  0.7859  (from  Table  3.3.1) 
are  achieved.  The  optimum  u*  differs  greatly  for  the  maintenance- 
dependent  A  From  the  u*  's  for  constant  A's.  This  indicates  that 
if  the  failure  rate  parameter  is  not  known  exactly,  the  characteristics 
of  the  system  changes  greatly  when  the  A 1 s  are  actually  maintenance- 
dependent  rather  than  constant. 


-‘i 


Fig.  4.4.9  p{a  J 1-  f  \  versus  u  for  various  failure  rate  parameters, 
c  * 


(4)  VARIATION  OF  REPAIR  RATE  PARAMETER  (p) 


Fig.  4.4.10  shows  the  variation  of  the  repair  rate  parameter  : 
p  »  1.0,  1.5,  and  2.0  (columns  (2),  (7),  and  (8)  in  Table  4.4.2). 

The  optimum  operative-cycle  time  u*  does  not  seem  to  change  ,  probably 
because  the  system  failure  rate  is  statistically  independent  of  the 
repair  time  T^  .  The  smaller  the  jx  the  longer  the  mean  repair  time, 
hence  the  smaller  the  availability.  Thus  the  probability  curves  in 
Fig.  4.4.10  are  lower  for  smaller  fi  . 


(5)  VARIATION  OF  SYSTEM  DESIGN  PARAMETER  (£) 

P  { 0.9J  and  p|ac^0.35J  versus  u  are  plotted  in  Fig.  4.4.11 

(columns  (2)  and  (9)  in  Table  4.4.2).  The  smaller  the  £  the  larger 

the  A  ,  hence  the  smaller  the  probability*  If  A  is  large,  a  long 
c  c 

optimum  operative-cycle  time  u*  has  to  be  used.  Hence  u*  is  larger 
when  £»  0.1  than  when  £  a  0.15  . 


a  •  2 

^  “  o.i 

F  m  1 
d  «  0.01 

g  »  0.5  (fixed  reduction) 


EXAMPLE  4.4.3 


The  age  replacement  policy  in  §  4.2  will  be  illustrated  in  this  example 
for  the  following  case  : 

d*  2,  £»  0.1,  0.1,  »  1,  d  ■  0.01,  g  -.©.S  (fixed  reduction) 

The  values  of  P^A^l-cJ  ,  tR  (the  age  replacement  time),  and  pj  A^  Jl-£ 
versus  the  operative— cycle  time  u  are  computed  the  PASCAL  computer  program. 
.  t 

The  graph  of  P  j  A^  £  l-£  J  versus  t^  for  some  values  of  u  are  plotted 

in  Fig.  4.4.12.  Table  4.4.3  depicts  the  maximum  P  ^  £  l-£ |  . 

The  age  replacement  policy  is  justified  by  comparing  Fig.  4.4.2  and 
Fig.  4.4.12. 


o»-  2 

Fixed  reduction,  g-0.5 

£=»  0.1,  a  0.1,  u  «1,  d  >  0.01 

p{ac>A-£| 

< 

maximum 

0.5 

0.5899 

25.49 

0.5904 

1 

0.G2B7 

28.27 

0.6287 

2 

0.G368 

32.15 

0.6368 

5 

0.G162 

35.06 

0.6162 

10 

0.5896 

40.03 

0.5896 

Table  4.4.3  Some  optimum  age  replacement  times. 


versus  t_  for  various 


1 


U4 


The  existence  of  on  optimum  age  replacement  time  (Fig.  4.4.12  and 

Table  4^4.3)  may  be  explained  by  the  following  argument.  Suppose  the  system 

does  not  fail  at  the  end  of  the  0 — I-i  cycle.  Tho  probability  of  system 

“tell 

failure,  in  the  (K+l)  0-M  cyclo  is  higher  than  in  the  previous  0— M  cycles 

J.L 

because  of  aging.  Hence  to  replace  the  system  at  the  end  of  the  It  1  0-cycle 
is  more  preferable  than  to  risk  over  the  next  0-M  cycle.  Thus  we  expect  an 

optimum  t*  for  smaller  operative-cycle  time  u  .  When  u  is  large,  the 
system  may  fail  during  the  first  few  0-M  cycles.  In  this  case,  t*  should 
be  chosen  as  the  first  time  that  the  probability  of  the  stochastic 
availability  achieves  certain  acceptance  level. (see  (4.2.14)}. 


NOTE  : 


The  time  scale  in  the  examples  is  not  necessary  the  real  time 
unit.  Interpretation  is  required.  For  example,  we  have 


0.1 


and 


Then 


l  u 


u 


0.2 

0.2/  X 


& 


(0.2256)- 

2  X 

( 0 . 2256 ).(  MTBF ) 


For  a  real  system  with  MTBF  ■  1000  hours,  say,  then  the 
optimum  operative-cycle  time  in  this  case  is 

u*  m  225.6  hours. 


Other  parameters  and  times  are  interpreted  similarly 


CHAPTER  5 


It  is  therefore  reasonable  to 

define  the  generalized  stochastic  availability 

asociated  with  costs  as 

T. 

A(°d«cf) 

3 

U 

(5.1.2) 

TU  *  °<1T0  *  cfTF 

m 

1  J  C(cd,c^J 

(5.1.3) 

Note  that 

(  0<C<1  , 

o<6<i) 

P 

f A(cd,cf)  *  1 

-4 

3  p  1  c^°d,cf^  ^ 1  +  rrr } 

(5.1.4) 

or. 

■  ?  ^ 1  ~  rfr  l 

Q 

i'c(cd,cf)  1 

♦jj 

(5.1.5) 

e- 

o 

1 

therefore  define  various  types  of  generalized  stochastic 

availabilities,  or  equivalently,  stochastic  cost  functions  as  in  §3.1 
and  §4.2. : 


(1)  Generalized  Stochastic  Cycle  Availability 
(Stochastic  Cycle  Cost  Function) 


VW 


operative 


operative 


c  T 

d  maintenance 


*  cfTF 


(5.1.0) 


where 


T 

operative 

T 

maintenance 


operative  uptime  for  a  renewal  cycle 
maintenance  downtime  for  a  renewal  cycle 


T. 


repair  time  for  a  renewal  cycle 


(Stochastic  Cost  Function  at  a  Finite  Time  in  the  n  Q-M  Cycle) 


VvV 


C  (c  ,c  J 
n  a  I 


n-i  A  n-x 

>  u.  +  rain(t,T  .  )  c  d .  +  c_T_ 

jz I  J  utn  ^  f  F 


,  0<  t«t-t  ,  <  u 
n-lx  n 


(5.1.10) 
,  0<t=t^!Tn^dn 


ZZUj  +  cd<  ZZdj  +  min(t,T  ))+  c  Tp 
j-1  J  j«l  J 


(6)  Generalized  Stochastic  Availability  for  an  Age  Replacement  Time 
(Stochastic  Cost  Function  for  an  Age  Replacement  Time) 


AAH(cd,Cf}  3  t„ 


CAR(°d’Cf) 


>  u.  +  T  , 

fe  J  U‘n 


- i - ^ - 7 -  *  T*  ftn  i'tn) 

n-1  n-1  n-l  n 

>  u.  +  T  .  +  c  .(>  d.  +  T, .  )  +  c„T  .  5  ,, 

*rrr  j  u  I  n  d  r—r  j  d  |  n  f  F  n=l ,  2 , . . . ,  K 

J*1  J"! 


>  u .  +  c  ,  >  d .  +  c  'r„ 

T- 1  J  di^  J  fr 


(5,1.11) 


Mote  that  when  c^  ■  cp  *  1  ,  the  above  definitions  reduce  to  those  defined 
in  Chapter  3  and  §4.2  . 


T 
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§5.2  PRODADILITY  DISTRIBUTIONS  OF  GENERALIZED  STOCHASTIC  AVAILABILITY 


•'  .'j 


The  probability  distributions  of  various  generalized  stochastic 
availabilities  or  stochastic  cost  functions  can  be  derived  in  the  same  way 
as  in  Chapter  3  and  |  4.2  .  We  shall  depict  Fig.  5.2.1  and  Fig.  5.2.2 
as  parallel  to  Fig.  3.2.1  and  Fig.  3.2.2  respectively  for  the  domains  of 
definitions  of  the  corresponding  variables,  ans  we  shall  summarize  the 
expressions  in  the  following  : 


0<l<  1 

U  =»  (  n  f  (z)  dz 

n  it  ,  n 
n-1 

°» -  Cn  f„(z)  - 


n— 1 

*  {*”  f„‘z)  - 


c.  +  f/(l-«)  n-1 

Z  -  -4-3 - DjI  o* 

f  j-1  J 


d  cd  +  £/(l-£)  rwl 

z  +  - ~ 

f  f  J 


n— j. 

z: )  dz 

j-i  J 


r  cd  +  1/  ( l-£ }  rwl 

tn-l  3  max  [  Vl  *  H  (1-0 


•c  f  -d 

un  3  max  I  min  (t,t  )  ,  — 


t  a  max 
n 


c.  +  !/(!-«)  n-1 


n— x 

H‘jl 

jal  J  ' 

n-1  -j 

*/(!-«)  dj  J 

{r  cd  +  n  -  -» 

min  [  min(t,tn)  ,  ,/^y  II  uj  J  ’ 


^n(cd,cf}  3  ttn^d’V] 
^n(cd,cf)  5  *flnlcd,cf,| 


tat 


tat 


t 

l 


piAc(cd‘cf^W  \  “  Un^»(Cd,Cf)  *  Dn^n(cd*0f) 


U  +  D 
n  n 


'IV^d-V)1'4}  -  [  UX(=d’V  *  “A(c d>V  J  1  (Un 

p{*0  (cd-ef>  *  ■  ^‘VV  /  <V  V 

n 

P  I  A«  <  Vcf>  ?  l-  61  -  4<  W  ‘  lUn  *  V 

n 

|  *hVf>  '  lun  *  °„>  • 
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CHAPTER  G 


CONCLUSIONS 


A  general  repairable  system  model  with  a  periodic  maintenance  schedule, 
and  age  and  maintenance  dependent  failure  rate  has  been  studied.  This  work 
examines  two  special  types  of  failure  reduction,  various  concepts  of 
stochastic  availabilities  and  their  probability  distributions,  optimum 
system  design  using  stochastic  availability,  age  replacement  policy,  failure 
rate  characteristics  and  parameter  variations,  computational  aspects,  os  well 
as  geneneralized  stochastic  availabilities  and  stochastic  cost  functions. 

For  a  periodic  operative-maintenance  schedule  with  failure  reduction, 

we  have  observed  the  existence  of  a  unique  optimum  operative-cycle  time  u* 

(Fig.  4.4.2  and  4.4.6).  An  analytical  proof  of  the  existence  and  uniqueness 

of  u*  is  difficult  (Eq. (3.2.18) ) .  However,  we  have  studied  the  asymptotic 

behavior  of  the  p{  A^l-ij  versus  u  curve  (§  3.3  and  Fig.  4.4.2).  We 

conclude  that  the  curve  with  failure  reduction  is  above  the  asymptotic  value 

P  and  hence  above  the  curve  with  no  reduction.  The  reason  is  that  the 
CO 

lifetime  probability  density  function  with  failure  reduction  is  asymptotically 

greater  than  the  p.d.f.  with  no  reduction  (Eq. (2.5.16)  ).  Thus  the 

P  f  A  y,  1-  f  \  crosses  the  P  asymptotically  from  above .  Hncce  the 
i  c '  ]  oo 

existence  of  an  optimum. (not  necessary  unique)  u  is  concluded.  As  a 
further  extension,  the  same  problem  with  non-periodic  operative— Maintenance 
policy  may  be  studied  to  provo  the  existence  and  possibly,  uniqueness  of  a 
set  of  optimum  operative-cycle  times. 

Some  optimum  system  design  criteria  have  been  proposed  and  illustrated 
by  examples  in  Chapter  4  to  obtain  an  optimum  oprative-cycle  time  and  an  age 
replacement  time.  'With  an  appropriate  interpretation  of  the  time  and 
parameter  scales (  as  noted  at  the  end  of  §  4.4  ) ,  we  can  use  the  results  and 
curves  in  §  4.4  in  system  design  to  achieve  an  optimum  schedule  and  to  study 
system  characteristics  against  parameter  variations.  Since  the  analysis 
involves  complicated  expressions,  the  PASCAL  computer  program  has  found  to 
be  useful  in  generating  results  as  well  as  insights  into  the  problem. 
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Some  parameter  variations  of  the  repairable  system  have  been  explored  in 
Example  4.4.2.  For  a  real  system  the  statistics  of  the  failure  rate  are 
usually  unknown  or  can  be  estimated  approximately  only.  The  sensitivity  of 
optimum  design  depends  on  parameter  variations.  In  particular,  the  optimum 
operative-cycle  time  u*  changes  greatly  for  a  maintenance-dependent  failure 
rate  parameter  A  than  constant  A's.  This  affect  is  important  in  system 
design  because  this  case  represents  a  system  with  a  rapid  varying  failure  rate 
and  it  is  therefore  different  from  the  constant  failure  rate  parameter  case. 
Other  parameter  variations  can  be  studied  in  a  similar  way.  For  example,  the 
repair  rate  u  needs  not  be  fixed,  or  the  repair  time  may  not  be  exponentially 
distributed,  or  the  repair  time  may  depend  on  the  degree  of  system  failure. 

As  possible  extensions  and  generalizations,  the  following  topics  are 
proposed  : 

(1)  Other  types  of  failure  reduction  criteria  and  nonlinear  failure  rates 
(not  necessarily  of  V/eibull  type)  can  be  defined  in  a  similar  way. 

(2)  More  general  repair  time  distribution  may  be  used.  For  repair  time 
depends  on  system  failure,  the  analysis  involves  the  evaluation  of  double 
integrals  instead  of  single  integrals.  This  can  be  done  both  theoretically 
and  numerically. 

(3)  Non-periodic  operative-maintenance  policy  may  be  useful  to  obtain  an 
optimum  set  of  operative-cycle  times  by  probabilistic  maximization  of  the 
stochastic  cycle  availability  when  the  system  parameters  change  rapidly 
with  age  and  maintenance. 

(4)  Age  replacement  policy  with  non-periodic  operative-maintenance  schedule 
may  be  treated  similarly. 

(5)  An  appropriate  interpretation  of  system  costs  will  apply  the  use  of 
generalized  stochastic  availabilities  or  stochastic  cost  functions  in 
system  design. 

(6)  The  present  PASCAL  computer  program  can  be  modified  to  satisfy  all  the 
above  extensions  and  potential  applications,  because  the  program  can  bo 
implemented-,  with  slight  changes,  in  all  computers  with  the  standard 
U.C.S.O.  PASCAL  language. 


APPENDIX  1 


WEIDULL  DISTRIBUTION 


The  probability  density  function  f(t)  of  a  Vioibull  distribution  ia 
given  by  (t'eibull  (1051)) 

f(t)  »  ^Oft*-1  expC-At*)  ,  t  ^0  (Al.I) 

whore 

A  *  scale  parameter 
U  a  shape  parameter 

The  shape  of  f(t)  is  sketched  in  Fig.  Al.I. 

When  or*  i,  f(t)  is  an  exponential  p.d.f. 

When  Of  »  2,  f(t)  is  a  Rayleigh  p.d.f. 


by 


The  distribution  function  F(t)  of  the  Weibull  distribution  is  given 


F(t)  -  )  f(t)dt 

*  0 

f  1  -  exp  (-  AV)  ,  t  £  0 
■  to  ,  t<0 


(A1.2) 


f(t) 


Fig.  Al.I  Weibull  probability  density  function*. 
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APPENDIX  2 


PASCAL  COMPUTER  PROGRAM 


Tho  PASCAL  computer  program  described  in  §  4.3  to  evaluate  various 
probability  distributions  of  stochastic  availabilities  is  listed  below. 

In  the  APPLE  II  PLUS  microcomputer  system  the  program  is  called  tho 
"SYSTEM. URIC. TEXT"  file  and  is  stored  in  a  5/j"  x  5J£"  floppy  disk  named  the 
'APPLE#: '  disk.  To  start  the  system,  a  disk  named  'APPLE3: '  is  inserted 
into  the  disk  drive  then  followed  by  the  'APPLE#:1  disk.  One  simply  types 
"R"  to  run  the  program.  The  program  can  be  implemented  in  other  computer 
systems  with  minor  modifications  because  the  PASCAL  system  is  the  standard 
U.C.S.D.  PASCAL  language. 
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.  KUCK.'iil  UlAi-WAU  IMP!)  I  ,  UU  i  PU  !'  )  i 

{  ****^.i,t**»**W*;li*«:k:x:)t.K.'fc***:X*********#******.l:***^*******^:(i«**.'tiA.ii4..k  ) 

<**  PROBABILITIES  uF  STOCHASTIC  AVAILABILITY  OF  A  REPAIRABLE  ** ) 
<**  SYSTEM  WITH  AGE  AND  MAINTENANCE  DEPENDENT  ** > 

<**  FAILURE  RATES  AND  OPTIONAL  AGE  REPLACEMENT  POLICY  ** > 

(**  APRIL  14 t  1982  ** > 

( %%****%%xxx*****%%T.w***%*x***%*%**X*%*%**%t*%*******%*X***%**%**  ) 
( *$S+*  ) 

USES  TRANSCEND? 

CONST  CONV  »  1.0E-7? 

51  ='  N  7? 

52  a7  PCAC>1-E37  ? 

53  a'  ECN<T)3  7? 

54  a'  ART  '  t 

55  a'  PC AAR 3  'I 

56  ay  PCA03  7? 

57  a'  PC  AM  3  7  ? 

58  a7  PC  AN  3  7  ? 

59  ay  PCN<T)»N37? 

TYPEO  a'<*#  NO  FAILURE  REDUCTION  **>7  ? 

TYPE1  a'<#*  FIXED  REDUCTION  **>7  ? 

TYPE2  =' <**  PROPORTIONAL  REDUCTION  **>' ? 

PA  =7  PC  SCA>7  * 

TIME=7  OPERATIVE-CYCLE  TIME  =  '? 


VAR  DN»UN»UD» ID » IU»T1 »T2»T3»T  *TO*UPO»R»RT  »RR»RF»RU»FO* 

TUP  r TDOUNt TUT»TDT  * F»FU» FD»DTUP*PRQB»PR»DP» RD»RO»G> 
EN * MUT, MLT * E »E1 * EE* DELTA* FIX, SCALE »PARAM*UPT, ART* 
AO * AM , AN , PN , AAR , MU » LAMDA * LAMAF ,LINF * LZERO t  RE AL ? 
LUB* N» ALPHA, ALFAl,POINT»TERM*KIND»K»L» 

SELECT , PRINT ALL*  INTEGER? 

EXIT.*  BOOLEAN? 

CHOICE* date;  STRING? 
filer;  interactive? 


PROCEDURE  QUTPUTCONTROL  ? 

BEGIN 

WRITELN< 7 OUTPUT  MEDIUM  l'  )? 

WRITELN< 7  1  =  CONSOLE7  )? 

WRITELN( 7  2  =  PRINTER7)? 

WRITE( 7  WHICH  NUMBER?  7  )?  READLN(  SELECT  )? 

CASE  SELECT  OF 

i;  choice;=7comsole;7 ? 

2;  choice;=7printer:7 ? 

END? 

IF  <  SELECT=2  )  THEN 
BEGIN 

URITE<  7  PRINT  ALL  ?  <YES=1*  N0=0  )  7  )?READLN< PRINTALL  ) 
END? 

END?  <*  OUTPUT  CONTROL  *> 


FUNCTION  POWER<X;  REAL?  YI  INTEGER);  REAL? 
VAR  Z;  REAL? 

BEGIN 

Z*.»1.0? 


UHiLh  Y  .  0  DO 
BEGIN 

WHILE  NOT  ODD(Y)  DO 
BEGIN 

Yt=Y  DIV  2i 

x:*x*x 

end; 

y:=y-i; 

z:=x*z 

end; 

power:=z 

end;  <*  INTEGRAL  POWER  *) 


FUNCTION  MAX<A:  REAL;  BJ  REAL):  REAL; 
BEGIN 

MAX :  =0  *  5*<  A+B+ABS<  A+B  )  ) 

END;  <*  MAXIMUM  *) 


FUNCTION  MINCA:  REAL;  BJ  REAL):  REAL; 
BEGIN 

min:=o.5*<  a+b-abs<  a-b  ) ) 
end;  <*  minimum  *) 


FUNCTION  EXPF<T:  REAL):  REAL; 

var  x:  real; 

BEGIN 

IF(T<-60.0)THEN  X.‘=0.0 

ELSE  X:=EXP<T); 

expf:=x 

end;  <*  EXPONENTIAL  FUNCTION  WITH  UNDERFLOW  SUPPRESSION  *) 


FUNCTION  REPAIR(T:  REAL):  REAL; 

BEGIN 

repair:=i.o-expf< -mu*t  )  i 

end;  <*  EXPONENTIAL  REPAIR  TIME  CDF  *) 


PROCEDURE  DISTRIBUTION TfTOfRRfDEL:  REALGAR  R»RF»F:  REAL); 

BEGIN 

R  : =PARAM*POWER<  T , ALF A1  )-DELTA ; 

RF:=EXPF< -LAMAF*< POUER<  T» ALPHA  )-POUER<  TO » ALPHA ) )+DELTA*<  T-TO  )  ) 
F  :=R*RF*RR 

END;  <*  AGE  &  MAINTENANCE  DEPENDENT  WEIBULL  FAILURE  RATE  *) 


PROCEDURE  SIMPSON< Tl»Il»T2f I3rRRfDELTA»EEfTTU»TTD»HUfHD:  REAL  J 

npu»npd:  integer;  var  iu»id:  read; 

VAR  SUO  f  SU1 >  SU2 » SU4 f  SDO » SD1 » SD2  *  SD4 » PDF  r  SF » 

X»Y»U1»U2»D1»D2:  real; 

j»nui»ndi:  integer; 

BEGIN 

iu:=o.o; 

IF  HU>0  *0  THEN 
BEGIN 

U1  :=MAX(0.0»EE*I1-TTD  )t 
DISTRIBUTION  1 1 » T 1 » RR  f  BELT  A  f  R  f  SF  f  PDF  )! 


suo:-puF*Rt:i'»iiK(  ui >; 

tJ2  :=EE*T2-TTD» 

DISTRIBUTION  T2  >  T 1 » RR » BELT  A  t  R » SF » PDF  )» 
SUi:=PDF*REPAIR<  U2)» 

su2:=o,o; 

SU4i=0.0» 

NUi:=NPU-l» 

FOR  JJ=1  TO  NU1  BO 
IF  QDD<  J  )  THEN 
BEGIN 

X  :=ii+j*hu; 

Y  :=EE*X-TTD» 

DISTRIBUTION X»Tl,RR»DELTA»RfSF, PDF) 
SU4i=SU4+PDF*REPAIR<  Y  ) 

END 

ELSE 

BEGIN 

X  :=ii+j*hu; 

Y  :=EE*X-TTDf 

DISTRIBUTION  X f T1 »RR > DELTA »R>SF»PDF > 
SU2J=SU2+PDF*REPAIR<  Y  ) 

end; 

IU :  =«HU*<  SU1+SU0+2 . 04SU2+4 . 0*SU4  )/3 , 0  i 

end; 

id:=o.o; 

IF  HD>0.0  THEN 
BEGIN 

B1  :*-T2+ttu; 

DISTRIBUTION  T2» Tit RR»DELTA*R»SF»PDF  )»* 
SDO:=PDF*REPAIR<DI ); 

D2  :=-I3+ttu; 

DISTRIBUTION  13 > T1 »RR» DELTA »R>SF> PDF  )f 

SDi:=PDF»REPAIR< D2H 

SD2J=0.0; 

SD4J=0.0; 

ndi:=npd~i; 

FOR  j:=l  TO  ND1  DO 
IF  ODD(J)  THEN 
BEGIN 

x  :=t2+j*hd; 
y  ;=-x+ttu; 

DISTRI BUTI ON  X » T1 r RR » DELTA »R»SF» PDF  ) 
SD4J=SD4+PDF*REPAIR<  Y ) 

END 

ELSE 

BEGIN 

x  :=t2+j*hd; 

Y  i=-X+TTU; 

DISTRIBUTION  X»T1»RR»DELTA»R»SF»PDF  ) 
SD2J=SD2+PDF*REPAIR<  Y ) 

end; 

IDl=HD#(  SD1+SD0+2.0*SD2+4.0*SD4  )/3.0» 

end; 

end;  <*  SIMPSON'S  RULE  *> 


PROCEDURE  INTEGRATION TlfT2»T3rRRf DELTA?  REAL; 

VAR  UPTflUflDJ  REAL); 

VAR  HD»HU»I1»I3»TTU»TTD»UTT:  real; 

np»npu»npd:  integer; 


BEGIN 

NPU:=2*R0UND< PQINT/2  )f 

npd:=2*round<  point/4  )i 

II  :»HIN<  T2fT1+MAX( 0 *0»-TUT+TDT/EE  ) )» 

utt:=tut+t2-ti; 

13  J=HAX<  0*OfMIN< T3fT2+EE*UTT-TDT  ) )» 

TTU:*UTT/E1 f 

ttd:*tdt/ei; 

HU  :»HAX(0*0»<T2-I1  )/NPU)» 

HD  J=MAX<0.0»<  I3-T2)/NPD)J 

SIMPSON T1  f I1fT2f I3fRRfDELTAfEEfTTUfTTDfHUfHDfNPUfNPDf IU » ID  ) 

upt:=utt; 

IU  J»HAX<  O.OfIU  )f 
ID  :=MAX<O.OfID) 

END;  <*  INTEGRATION  *) 


PROCEDURE  WEIBULL<Ti:  REAL ;  VAR  ROJ  REAL); 
BEGIN 

lamda:=linf-<  linf-lzerq  >*exp< -< n-i >/scale  >; 

ALFA1 i=ALPHA-l i 

LAMA F : =POWER< LAMDA f ALPHA  ) ! 

PARAM : =ALPHA*LAMAF ; 

RO  :=PARAM*POUER<  T1 fALFAI  ) 

END;  <*  UEIBULL  PARAMETER  *) 


PROCEDURE  PROBABILITY; 

BEGIN 

pr  :=prob; 

DISTRIBUTION T2fT1 fRRfDELTAfRfRUfFU  )» 
DISTRIBUTION  T3 f TIfRRf DELTA fRT fRDfFD  >; 
un  :=max<o.OfRR*<  i.o-RU)>; 

DN  :=MAX<0.0»RR*<RU-RD))f 
ud  j=dn+un; 
pn  :=pn-ud; 
fo  :=fo+un; 

INTEGRATION  T1fT2fT3fRRf DELTA f UPT f IU  f IB  ) ; 

IF  <UD>0.0>  THEN 
BEGIN 

ao  :=iu/ud; 

AM  5=ID/UD; 

AN  :=UN*AO+DN*AMf 

pros  :=prob+an; 

AN  := AN/UN 

end; 

art  :*upt+tdt; 

AAR  :=PR+IU+PN*REPAIR<  UPT/E1-ART  )i 

en  :=en+n*ud; 

IF  <  PRINTALL=*1  )  AND  <  SELECTS  )  THEN 

WRITELN  FILER f N  J3fPR0B » 10 1 6 f EN  1 7 1 2  f  ART ♦ 8 J3f 

aar:9:6fAo:io:6fAm:io:7fAn:io:7fUd:i2); 
WRITELN NJ3fPR0B : 10 J6fEN :7J2f ART  58:3 F AAR :?J6f 

ao : io:6f am: io :7f an: io :7 fUd: i2>; 
term  :=n; 

TUT  :=TUT+T2-T1f 
TDT  :»TDT+T3-T2; 

rr  :=rr*rd; 

ti  :*T3; 

T2  :»ti+tup; 


1.4  ;  -  iLMTDOUN  > 

n  :=n+i; 

WEIBULL*  T1 , RO  ) » 

IF  <  KIND=1  )  THEN 

DELTA  :=MAX*0.0,R0-RT+FIX >; 

IF  *  KIND-2)  THEN 

DELTA  J  -MAX*  0 » 0 , RO-< 1 » O-G  )*RT  ) ; 

DP  J=ABS*PROB-PR); 

IF  *DP<CONV)  AND  <PRQB<1.0)  AND  <N>1>  THEN  EXITt=TRUE 
END?  **  TOTAL  PROBABILITY  *) 


PROCEDURE  INPUT? 
BEGIN 

uriteln; 


WRITE* 7  **** 
WRITE( 7 
WRITE* 7 
WRITE*  7  LAMDA 


7  LAMDA 

:-io.o; 


WRITE* 

SCALE 
WRITE( 7 
WRITE( 7 
WRITELN* 
WRITELN* 
URITELN< 
WRITELN* 
WRITE*  7 
IF  <KIND=0) 
ELSE 
BEGIN 
WRITE* 


DATE  J 

E 

ALPHA 

INITIAL 

FINAL 


>;READLN(  DATE); 

7  );readln<  e>; 

7  >;READLN(  ALPHA); 

7  );readln<  lzero); 
7  );readln(linf); 


mu  = 


>;readln(  mu); 

MAINTENANCE-CYCLE  TIME  =  7  ) ; READLN( TDOWN  ) ; 
CRITERION  OF  FAILURE  REDUCTION  J  7  ); 

TYPE  0  *  NO  FAILURE' REDUCTION  '  )i 
TYPE  1  =  FIXED  REDUCTION  7  ); 

TYPE  2  a  PROPORTIONAL  REDUCTION  7  >; 

type  =  7  );readln(kind>; 

THEN  G:=0.0 


REDUCTION  FACTOR  =  7  );READLN(G); 


end; 

WRITE*  7  INITIAL  O-CYCLE  TIME  = 
O-CYCLE  TIME  INCREMENT  = 


WRITE* 

WRITE* 'NUMBER  OF  O-CYCLE  TIME  = 
URITE* 7  SUMMATION  TERMS  = 
WRITE* 7  INTEGRATION  POINTS  = 

END;  (*  INPUT  DATA  *) 


);readln<  upo); 
);readln<  dtup); 
);readln<l>; 
);readln<  lub); 
);readln<  POINT) 


PROCEDURE  DATA; 

BEGIN 

input; 

lamda:=lzero; 
ei  :=i*o-e; 
ee  :=e/ei; 
rewrite*  filer , choice >; 

WRITELN<  FILER); 

URITELN* FILER. 7  DATE  J  7  ,DATE  );WRITELN<  FIL  ER  ); 

UJRITELN*  FILER, 7  STOCHASTIC  AVAILABILITY  **7  ); 

IF  ( KIND**)  )  THEN  WRITELN*  FILER, TYPE*) ) 

ELSE 

BEGIN 

IF  (KIND-1)  THEN  URITELN* FILER, TYPE1  ); 

IF  <  KIND=2  )  THEN  URITELNt FILER, TYPE2 >; 

WRITELN*  FILER, 7  REDUCTION  FACTOR  =  7,Gi7:5> 

end; 

URITELN* FILER,7  E  =  7,E:5{2>; 

URITELN* FILER,7  APLHA  =  7,ALPHAJ2>; 


J.w 


LAMDA 


t  L  Aii£>A  1714  ) 


ishiTELM:  FILER i' 

EL 'Jc 
BEGIN 

URITELM FILER* 7 LAMDA  INITIAL  = 
WRITELN( FILER* 7  LAMDA  FINAL  = 
END  i 


7  *LZERQ»7»4 >} 
7  »LINFJ7 *4  ) 


URITELM  FILER* 7  MU  =  7»MUJ7:4)» 


URITELM  FILER* 7  MAINTENANCE-CYCLE  TIME  =  7 »TD0UN:7i5  )> 

WRI »ELN<  FILER* 7  INTEGRATION  =  7 .POINT*7  POINTS7) 


END»<*  PRINT  DATA  4) 


iivi 


PROCEDURE  HEADING r 
BEGIN 

IF  <PRINTALL=1>  AND  < SELECT=2 >  THEN 
BEGIN 

WRITELNI FILER )J 

UR I TELNt FILER » T IME » TUP  J 8  J 5  ) ; 

URITELNC FILER . S 1 #  S2 f  S3 . S4 * S5 * S6 * S7 *  S8 » S9  ) 
END? 

WRITELN<  )i 

URITELM  7  .  ♦  ♦  ♦  .7  »TIME*TUPt8*5  )» 

URITELM  51 *S2*S3*S4*S5*S6*S7*S8*S9  ); 

URITELM  ) 

END;  <*  PRINT  HEADING  #> 


PROCEDURE  INITIAL; 
BEGIN 


EXIT 

: -false; 

ALFA1 

:=alpha-i; 

LAMAF 

J  =POUER<  LZERQ , ALPHA  > ; 

PARAM 

:=alpha*lamaf; 

FIX 

:=G*PARAM*<  TUP+TDOWN  >; 

DELTA 

:=o.o; 

N 

H 

II 

PROD 

:=o»o; 

PN 

:=i.o; 

RR 

;=i.o; 

FQ 

:=o,o; 

TDT 

:=o.o; 

TUT 

;=o,o; 

EN 

49k 

o 

♦ 

o 

II 

T1 

II 

o 

• 

O 

*•* 

T2 

:=tup; 

T3 

:=T2+TB0WN 

end;  <# 

INITIAL  conditions  *) 

BEGIN  <**  MAIN  PROGRAM  ** ) 

outputcqntrol; 

data; 

tup:=upo; 

FOR  K*=l  TO  L  DO 
BEGIN 

initial; 

heading; 

REPEAT  PROBABILITY  UNTIL  <  N>LUB  )  OR  EXIT! 
IF  N>LUB  THEN 


BEGIN  105 

WRITELN( FILER >» 

WRITELN< FILER,'**  SCA  NOT  CONVERGE  IN  ' »LliB» '  TERMS  **'  )i 
WRITELN<  FILER#PA»£1 J4t2»'  3  -  7 >PRQBJ9J6» ' »  ' >TIMEfTUP) 

END 

ELSE 

BEGIN 

hut:»tup*< en-o.s  >; 

MLT : =<  TUP+TDOWN  >*<  EN-0 .5  )i 

WRITELN<  FILER  )»URITELN< FILER » TIME »TUPJ8{ 5  )i 

WRITELN<  FIL£R»PA»E1 *4*2» '  1=  ' f?RQBl9l6f '  >  €  '  fTERM» '  J'  >; 

WRITELN<  FILER* ' MEAN  0-M  CYCLES  --  ' »EN  )» 

WRITELN( FILER* '  MEAN  UPTIME  =  ' *MUT >» 

WRITELN( FILER*'  MEAN  LIFETIME  =  '  »MLT  ) 

end; 

tup:=tup+dtup 

end; 

URITELNC  );WRITELN< ' *************  END  ******************'  ); 
WRITELNC FILER*' *************************************'  * 

'  *******************************************'  >; 

CLOSE<  FILER) 

END.  <**  STOCHASTIC  AVAILABILITY  **  > 
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